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Warm Dark Matter (WDM) research is progressing lastly, the subject is new and WDM essen- 
tially works, naturally reproducing the astronomical observations over all the scales, small (galactic) 
as well as large and cosmological scales (AWDM). Evidence that Cold Dark Matter (ACDM) and 
its proposed tailored cures do not work at small scales is staggering. The Chalonge Workshop 
'Warm Dark Matter in the Galaxies: Theoretical and Observational Progresses', was held at the 
Meudon Chateau of Observatoire de Paris on 8-10 June 2011. The Workshop approached DM in a 
fourfold coherent way: astronomical observations of DM structures (galaxy and cluster properties, 
haloes, rotation curves, density profiles, surface density and scaling laws), AWDM N-body simula- 
tions in agreement with observations, WDM theoretical astrophysics and cosmology (kinetic theory, 
Boltzmann-Vlasov evolution, halo model, improved perturbative approachs), WDM particle physics 
(sterile neutrinos) and its experimental search. Fedor Bezrukov, Pier-Stefano Corasaniti, Hector J. 
de Vega, Stefano Ettori, Frederic Hessmann, Ayuki Kamada, Marco Lombardi, Alexander Merle, 
Christian Moni Bidin, Angelo Nucciotti on behalf of the MARE collaboration, Sinziana Paduroiu, 
Henri Plana, Norma G. Sanchez, Patrick Valageas, Shun Zhou present here their highlights of the 
Workshop. Cored (non cusped) DM halos and WDM (keV scale mass) are clearly determined from 
theory and astronomical observations, they naturally produce the observed structures at all scales; 
keV sterile neutrinos are the most suitable candidates, they naturally appear in minimal extensions 
of the standard model of particle physics. AWDM simulations with keV particles remarkably repro- 
duce the observations, the small and large structures, sizes of local minivoids and velocity functions. 
The summary and conclusions by H. J. de Vega and N. G. Sanchez stress among other points the im- 
pressive evidence that DM particles have a mass in the keV scale and that those keV scale particles 
naturally produce the small scale structures observed in galaxies. Wimps (DM particles heavier than 
1 GeV) are strongly disfavoured combining theory with galaxy astronomical observations. keV scale 
sterile neutrinos are the most serious DM candidates and deserve dedicated experimental searchs 
and simulations. Astrophysical constraints including Lyman alpha bounds put the mass in the range 
1 < m < 13 keV. Predictions for EUCLID and PLANCK have been presented. Interestingly enough, 
MARE -and hopefully an adapted KATRIN- experiment could provide a sterile neutrino signal. It 
will be a a fantastic discovery to detect dark matter in a beta decay. There is a formidable WDM 
work to perform ahead of us, these highlights point some of the directions where it is worthwhile to 
put the effort. Photos of the Workshop are included. 
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I. PURPOSE OF THE WORKSHOP, CONTEXT AND INTRODUCTION 

This Workshop addresses for the first time the turning point in the research of Dark Matter represented by Warm 
Dark Matter (WDM) putting together astrophysical, cosmological and particle WDM, astronomical observations, 
theory and WDM numerical simulations which naturally reproduce the observations, as well as the experimental 
search for the WDM particle candidates (sterile neutrinos). 

Recently, AWDM emerged impressively over ACDM (A-Cold Dark Matter) whose small -galactic- scale problems 
are staggering. AWDM solves naturally the problems of ACDM and agrees with the observations at small as well as 
large and cosmological scales. 

This Workshop is the second of a new scries dedicated to Dark Matter. 

The first Workshop of this series in the Meudon Castle CIAS in June 2010 allowed to identify and understand 
the issues of the serious problems faced by Cold Dark Matter (CDM) to reproduce the galactic (and even cluster of 
galaxies) observations. 

The 2010 and 2011 Workshops served as well to verify and better understand the always growing amount of confusion 
in the CDM research, namely the increasing number and cyclic change of arguments, counter-arguments and ad-hoc 
mechanisms introduced in the CDM simulations over most of twenty years, in trying to deal with the CDM small scale 
crisis: Cusped profiles and overabundance of substructures. Too many satelites are predicted by CDM simulations 
while cored profiles and no such overabundant substructures are seen by astronomical observations. A host of ad-hoc 
mechanisms are proposed and advocated to cure the CDM problems. 'Baryon and supcrnovae feedbacks', non circular 
motions, triaxiality, mergers, 'cusps hidden in cores', 'strippings' are some of such mechanisms tailored or exagerated 
for the purpose of obtaining the desired result without having a well established physical basis. For example, the strong 
"baryon and supernovae feedback" introduced to transform the CDM cusps into cores in baryon+CDM simulations 
corresponds to a large star formation rate contradicting the observations. 

On the CDM particle physics side, the problems are no less critical: So far, all the dedicated experimental searches 
after most of twenty years to find the theoretically proposed CDM particle candidate (WIMP) have failed. The CDM 
indirect searches (invoking CDM annihilation) to explain cosmic ray positron excesses, are in crisis as well, as wimp 
annihilation models are plagued with growing tailoring or fine tuning, and in any case, such cosmic rays excesses are 
well explained and reproduced by natural astrophysical process and sources, as non-linear acceleration, shocks and 
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FIG. 1: Poster of the Workshop 



magnetic winds around massive explosion stars, quasars. The so-called and repeatedealy invoked 'wimp miracle' is 
nothing but been able to solve one equation with three unknowns (mass, decoupling temperature, and annhiliation 
cross section) within WIMP models theoretically motivated by SUSY model building twenty years ago (SUSY was 
very fashionable at that time and believed a popular motivation for many proposals). 

After more than twenty years -and as often in big-sized science-, CDM research has by now its own internal inertia: 
growing simulations involve large super-computers and large number of people working with, CDM particle wimp 
search involve large and longtime planned experiments, huge number of people, (and huge budgets); one should not 
be surprised in principle, if a fast strategic change would not yet operate in the CDM and wimp research, although 
they would progressively decline. 



FIG. 2: Photo of the Group 



In contrast to the CDM situation, the WDM research situation is progressing fast, the subject is new and WDM 
essentially works, naturally reproducing the observations over all the scales, small as well as large and cosniological 
scales (AWDM). 

This 2011 Workshop addressed the last progresses made in WARM DARK MATTER and the Universal and Non 
Universal properties of Galaxies. In the tradition of the Chalonge School, an effort of clarification and synthesis is 
made by combining in a conceptual framework, theory, analytical, observational and numerical simulation results 
which reproduce observations. 

The subject have been approached in a fourfold coherent way: 

(I) Conceptual context 

(II) Astronomical observations linked to the galaxy structural properties and to structure formation at different 
large and small (galactic) scales. 

(III) WDM Numerical simulations which reproduce observations at large and small (galactic) scales. 

(IV) WDM particle candidates, kcV sterile neutrinos: particle models and astrophysical constraints on them 

The Topics covered included: 

Recent progress in solving the Boltzmann-Vlasov equation to obtain the observed properties of galaxies (and clusters 
of galaxies). N-body numerical simulations with Warm Dark Matter; the surface density; scaling laws, universality 
and Larson laws; the phase-space density. Particle model independent analysis of astrophysical dark matter. The 
impact of the mass of the dark matter particle on the small scale structure formation. 



5 



The radial profiles and the Dark Matter distribution. Observed cored density profiles. The ever increasing prob- 
lems of ACDM at small scales. The keV scale Dark Matter (Warm Dark Matter): Observational and theoretical 
progresses. Perturbative approachs and the Halo model. Large and small scale structure formation in agreement with 
observations at large and small (galactic) scales. The new serious dark matter candidate: Sterile neutrinos at the keV 
scale: Particle physics models of sterile neutrinos, astrophysical constraints (Lyman alpha, Supernovae, weak lensing 
surveys), experimental searches of keV sterile neutrinos 

Fedor Bezrukov, Pier-Stefano Corasaniti, Hector J. de Vega, Stefano Ettori, Frederic Hessmann, Ayuki Kamada, 
Marco Lombardi, Alexander Merle, Christian Moni Bidin, Angelo Nucciotti on behalf of the MARE collaboration, 
Sinziana Paduroiu, Henri Plana, Norma G. Sanchez, Patrick Valageas, J. Shun Zhou present here their highlights of 
the Workshop. 

The Summary and Conclusions: 

Summary and conclusions are presented at the end by H. J. de Vega and N. G. Sanchez in three subsections: 

A. General view and clarifying remarks. 

B. Conclusions. 

C. The present context and future in the DM research. 

The conclusions stress among other points the impressive evidence that DM particles have a mass in the keV scale 
and that those keV scale particles naturally produce the small scale structures observed in galaxies. Wimps (DM 
particles heavier than 1 GcV) are strongly disfavoured combining theory with galaxy astronomical observations. keV 
scale sterile neutrinos are the most serious DM candidates and deserve dedicated experimental searchs and simulations. 
Astrophysical constraints including Lyman alpha bounds put the sterile neutrinos mass in the range 1 < m < 13 keV. 
Predictions for EUCLID and PLANCK have been presented. MARE -and hopefully an adapted KATRIN- experiment 
could provide a sterile neutrino signal. The experimental search for these serious DM candidates appears urgent and 
important: It will be a a fantastic discovery to detect dark matter in a beta decay. 

There is an encouraging and formidable WDM work to perform ahead us, these highlights point some of the 
directions where to put the effort. 

Sessions lasted for three full days in the beautiful Meudon campus of Observatoire de Paris, where CIAS 'Centre 
International d' Ateliers Scicntifiques' is located. All sessions took place in the historic Chateau building, (built in 
1706 by the great architect Jules-Hardouin Mansart in orders by King Louis XIV for his son the Grand Dauphin). 

The meeting was open to all scientists interested in the subject. All sessions were plenary followed by discussions. 
The format of the Meeting was intended to allow easy and fruitful mutual contact and communication. Large time 
was devoted to discussions. All informations about the meeting, are displayed at 

http://www.chalonge.obspm.fr/CiasJVreudon2011.html 

The presentations by the lecturers are available on line (in .pdf format) in 'Programme and Lecturers' in the above 
link, as well as the Album of photos of the Workshop: 

http:/ /www. chalonge.obspm.fr/Programme_CIAS2011. html 



http: / /www. chalonge.obspm.fr/albumCIAS2011/index.html 

We thank all again, both lecturers and participants, for having contributed so much to the great success of this 
Workshop and look forward to seeing you again in the next Workshop of this series. 

We thank the Observatoire de Paris and its Scientific Council, CNRS-GISP2I and UPMC supports, as well as the 
CIAS, secretariat and logistics assistance and all those who contributed so efficiently to the successful organization of 
this Workshop. 

With compliments and kind regards. 
Hector J de Vega, Norma G Sanchez 
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II. PROGRAMME AND LECTURERS 

• Fedor BEZRUKOV (Ludwig-Maximiliaiis-Universitat Miinchen, Germany) Light sterile neutrino dark matter 
in extensions of the Standard Model 

• Christopher J. CONSELICE (School for Physics and Astronomy, University of Nottingham, UK) Correla- 
tions between baryons and dark matter during galaxy formation 

• Pier Stefano CORASANITI (CNRS LUTH Observatoire de Paris, Meudon, France) Thcorical modeling of 
the halo mass function and halo bias: a path-integral approach to the Excursion Theory 

• Hector J. DE VEGA (CNRS LPTHE Univ de Paris VI, France) Warm Dark Matter from theory and galaxy 
observations 

• Stefano ETTORI (INAF, Osservatorio Astronomico Bologna, Italy) Mass profiles and concentrations in X-ray 
galaxy clusters 

• Frederic V. HESSMAN & Monika ZIEBART (Inst. F. Astrophysik, Georges- August-Univ. Gottingen, 
Gottingen Germany) The Bosma effect revisited - Correlations between the Interstellar Medium and Dark 
Matter in galaxies 

• Ayuki KAMADA (IPMU, Institute for the Physics and Mathematics of the Universe, Tokyo, Japan) keV-mass 
sterile neutrino dark matter and the structure of galactic haloes 

• Marco LOMBARDI (University of Milano, Italy) Larson' laws and the universality of molecular cloud struc- 
tures 

• Katarina MARKOVIC (Max-Planck-Ludwigs-Maximilians Universitt, Munich, Garching, Germany ) Con- 
straining Warm Dark Matter with future cosmic shear data 

• Alexander MERLE (Royal Institute of Technology KTH Stockholm, Sweden) Neutrino model building and 
keV sterile neutrino Dark Matter 

• Christian MONI BIDIN (Universidad de Concepcion, Departamento de Astromomia, Chile) The lack of 
evidence of a Dark Matter disk in the Milky Way 

• Angelo NUCCIOTTI (INFN Milano Biccoca, Italy) The MARE experiment in Milano Bicocca and its capa- 
bilities to measure the mass of light (active) and heavy (sterile) neutrinos 

• Sinziana PADUROIU, Observatoire de Geneve, Switzerland Numerical simulations with Warm Dark Matter: 
The effects of free streaming on Warm Dark Matter haloes of Galaxies 

• Henri Michel PLANA (Laboratorio de Astrofisica Teorica e Observacional, UESC Santa Cruz, Bahia, Brazil) 
Mass Distribution in Groups of Galaxies 

• Oleg RUCHAYSKIY (CERN-Teory, Geneva, Switzerland) Sterile neutrinos as dark matter candidates and 
the constraints on them from Lyman-alpha forest data 

• Paolo SALUCCI (SISSA-Astrophysics, Trieste, Italy) Structural properties of Galaxies and Cored Density 
Profiles from Galaxy observations 

• Norma G. SANCHEZ (CNRS LERMA Observatoire de Paris, Paris, France) Warm Dark Matter and Galaxy 
properties from primordial fluctuations and observations. 

• Joop SCHAYE (Leiden University, The Netherlands) Dark matter haloes, central AGN black holes, and the 
effect of baryons on the clustering of matter 

• Patrick VALAGEAS (Inst, de Physique Theorique, Orme de Merisiers, CEA-Saclay, Gif-sur-Yvette, France) 
Perturbation theories for large-scale structures 

• Shun ZHOU (Max-Planck-Heinsenberg Inst., Munich, Germany) Supernova bound on keV-mass sterile neu- 
trinos 
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III. HIGHLIGHTS BY THE LECTURERS 



A. Fedor Bezrukov 



Department fur Physik, Ludwig-Maximilians-Universitat Munchcn, Germany 
Institute for Nuclear Research of RAS, Moscow, Russia 

Light sterile neutrino dark matter in extensions of the Standard Model 

In this talk I would like to summarize the possibilities of having a light sterile neutrino DM in the extensions of the 
SM. Such a candidate is naturally very weakly interacting. At the same time, if sterile neutrino is light (of the order 
of several keV), its lifetime can exceed significantly the Universe lifetime without exceptional fine tuning. However, 
obtaining the proper observed abundance of the neutrino DM is quite constraining. The usual DM production 
mechanism is the tehrmal one, which corresponds to a particle that enters the thermal equilibrium and freezes out 
later. Then its abundance (number to entropy density ratio) is constant after freeze out (if no additional entropy 
is generated later), and depends on how relativistic the particle was at decoupling. If the particle freezes out while 
relativistic (what is the case for a keV scale particle), then its abundance is n/s = 135C,{3)/4^T^g^,f , with g^,f being the 
effective number of degrees of freedom at freeze out. Then the mass M of the particle is fixed from the requirement 
of its energy density being equal to the observed DM value 



1 /10.75\ / M 



^BM S \ (7*f / V 1 keV 

where S stands for the possible entropy generation factor. 
The possible situations are the following. 



X 100 , 




FIG. 3: Schematic evolution of the relic abundance in the Universe. Left plot depicts the standard CDM mechanism, where 
the abundance is controlled by the decoupling moment (several dashed curves). Right plot illustrates various non-thermal 
mechanisms. The dashed line is a thermal relic decoupled while being relativistic (hot thermal relic), leading to the over closure 
of the Universe. The blue decreasing line is the same hot thermal relic, but with the abundance diluted by rapid expansion of 
the Universe (entropy production), leading to correct DM abundance. The lowest (magenta) line depicts the evolution of the 
non-thermally produced particle with zero primordial abundance. 

Hot Dark Matter. If neutrino enters thermal equilibrium and there is no entropy dilution = 1, then the DM 
mass is M ^ 10 eV. This gives us HDM, which is excluded now from the structure formation observations. 

Entropy dilution. However, the constraint on the mass may be weakened if there is large subsequent entropy 
dilution S. This entropy production happens if there is some other particle with long lifetime, which first decouples 
while still relativistic and then decays when already non- relativistic [1]. Then the proper DM abundance is controlled 

by the properties of this long-lived particle through the entropy dilution factor S ~ 0.76 .'^^^ , where is an 

averaged number of d.o.f. during entropy generation. 

This situation can be realized if the sterile neutrinos are charged under some beyond the SM gauge group [2]. 
Note, that teh gauge scale should be below GUT, otehrwise the neutrinos do not enter the thermal equilibrium after 
reheating and the situation is similar to those of the i/MSM (described in the next paragraph) . In this setup one of 
the sterile neutrinos A'^i is light and is the DM particle, while other sterile neutrinos iV2.3 should dilute its abundance 
up to the correct amount by out-of-equilibrium decay. 
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FIG. 4: Bounds on the mass Mi and the mixing angle 6i of the DM sterile neutrino for (left to right) the model with entropy 
dilution, mtMSM, and production in the inflaton decays. 

The parameters of the DM sterile neutrino are constrained by two considerations. First, its lifetime in the radiative 
decays A^i — >■ 1^7 should be large enough not to contradict the X-ray observations. This constraints the mixing angle 
9i of A'^i and active neutrinos from above (see figure). At the same time the structure formation from the Lyman-a 
analysis constraints its mass from below Mi > 1.6 keV (its velocity distribution is that of the cooled thermal relic, so 
the bound is weaker, than in the other cases, analyzed in this note). 

All other constraints in this scenario correspond to the heavier sterile neutrinos and to the gauge sector. The 
right abundance of the sterile neutrino requires entropy dilution. To provide proper the entropy dilution, N2 should 

decouple while relativistic and has decay width F ~ 0.50 x lO^^^^gi^^^j^ ( ^ ) • same time, the heavy 

neutrino N2 should decay before BBN, which bounds its lifetime to be shorter than approximately 0.1 -^- 2 s. Then, 
the proper entropy can be generated only if its mass is larger than 

The entropy is effectively generated by out-of-equilibrium decay, if the particle decoupled while still relativis- 
tic. The bound on the decoupling temperature leads to the bound on the right-handed gauge boson mass 

A^>77.(r#vl)'^'(10-16)TeV. 

Non-thermal production within i/MSM. Another option to escape the overproduction of the light sterile 
neutrinos is to make them not enter the thermal equilibrium at all. This is the case for the simple extension of the SM 
by three right handed sterile neutrinos — i/MSM [3-5] . In this case the production of the sterile neutrinos proceed via 
oscillations from the active neutrinos. These oscillations are non-resonant if there are no lepton asymmetries during 
the DM generation, or resonantly enhanced if the lepton asymmetries are present. Then, the amount of produced 
DM is governed by the mixing angle Oi. However for nonzero lepton asymmetries, there is an allowed region [5,6] 
in 9i (see figure 2]), with the upper dashed line corresponding to the non-resonant production, and the lower shaded 
area corresponds to lepton asymmetries which are too large and contradict BBN. The usual X-ray bound also bounds 
9i. Together with the Lyman-a bound this forbids the non-resonant production mechanism, however there is an 
open region for the resonant mechanism. Note, that the Lyman-a bound depends strongly on the properties of the 
resonance, and can be as small as Mi > 2 keV for some values of lepton asymmetry, see [6] for details. 

Production in inflaton decays. Finally, the neutrinos may be produced by some other mechanism, like inflaton 
decay [7-9]. In this case the production is again not connected with 9i, so the only bounds present for the DM 
neutrino properties are the X-ray bound and the Lyman-a bound. It is interesting to note, that the mass of the DM 
sterile neutrino can be related to the mass of the inflaton, which is light in the model and can be searched for in 
B-meson decays [9]. 
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B. Pier Stefano Corasaniti & Ixandra Achitouv 

Laboratoire Univers et Theories (LUTh), 
UMR 8102 CNRS, Obscrvatoire de Paris, Univcrsite Paris Diderot, 
5 Place Jules Janssen, 92190 Mcudon, France 

Theoretical Modeling of the Halo Mass Function: A Path-Integral Approach to the Excursion Set 

Theory 

The computation of the halo mass function is of primary importance in cosmology. Halos are the building blocks 
of the cosmic structure formation and predicting their distribution over space and time is key to understanding the 
formation of galxies and the other visible structures in the universe. 

Because of the non-linear gravitational interactions which lead to the formation of halos, their properties have been 
mainly investigated using large volume numerical N-body simulations. 

The excursion set theory [1] provides a powerful mathematical framework to model the characteristics of the halo 
mass distribution. In the original formulation an analytical evaluation has been limited to few ideal cases, while 
requiring the use of Monte Carlo simulations to solve the stochastic model equations for more realistic scenarios. The 
recent formulation of the excursion set in terms of path-integrals [2] extends the analytical computation of the halo 
mass function to halo mass definitions which are consistent with those from observations, thus allowing for a direct 
comparison with numerical simulations as well as observational results. Furthermore, a stochastic modeling of the 
halo collapse conditions can be easily implemented in this formalism [3]. 

In this talk we review the derivation of the halo mass function in the case of a model of ellipsoidal collapse of halos 
[4,5]. We show that the inferred mass function can be written in terms of physically motivated quantities. Morevoer, 
the result of the theoretical computation is in excellent agreement with results from N-body simulation data. 
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keV scale dark matter from theory and observations and galsLxy properties from linear primordial 

fluctuations 

In the context of the standard Cosmological model the nature of dark matter (DM) is unknown. Only the DM 
gravitational effects are noticed and they are necessary to explain the present structure of the Universe. DM particle 
candidates are not present in the standard model (SM) of particle physics. Particle model independent theoretical 
analysis combined with astrophysical data from galaxy observations points towards a DM particle mass in the keV 
scale (keV = 1/511 electron mass) [1-4]. Many extensions of the SM can be envisaged to include a DM particle 
with mass in the keV scale and weakly enough coupled to the Standard Model particles to fulfill all particle physics 
experimental constraints. 

DM particles can decouple being ultrarelativistic (UR) at m or non-rclativistic T^i <^ m. They may decouple 

at or out of local thermal equilibrium (LTE). The DM distribution function: Fc;[pc] freezes out at decoupling becoming 
a function of the comoving momentum pc =. Pf(t) = pc/a{t) = is the physical momentum. Basic physical quantities 
can be expressed in terms of the distribution function as the velocity fluctuations, {V'^{t)) = {P'j{t)) /m? and the DM 
energy density PDAiit) where y ~ Pf{t) /Td{t) = Pc/Ta is the integration variable and g is the number of internal 
degrees of freedom of the DM particle; typically 1 < g < 4. 

Two basic quantities characterize DM: its particle mass m and the temperature Td at which DM decouples. 
Td is related by entropy conservation to the number of ultrarelativistic degrees of freedom gd at decoupling by 

Td = {2/gd)^ Tcmb , Tcmb = 0.2348 10^'^ eV. One therefore needs two constraints to determine the values of m and 
Td (or gd). 

One constraint is to reproduce the known cosmological DM density today. p£)A/(today) = 1.107 keV/cm^. 

Two independent further constraints are considered in refs. [1-4]. First, the phase-space density Q — p/cr^ [1-2] 
and second the surface acceleration of gravity (surface density) in DM dominated galaxies [3-4] . We therefore provide 
two quantitative ways to derive the value m and gd in refs. [1-4]. 

The phase-space density Q is invariant under the cosmological expansion and can only decrease under self-gravity 
interactions (gravitational clustering). The value of Q today follows observing dwarf spheroidal satellite galaxies of 
the Milky Way (dSphs): Q today ~ (0.18 keV)"* (Gilmore et al. 07 and 08). We compute explicitly Qprim (in the 
primordial universe) and it turns to be proportional to [1-4]. 

During structure formation Q decreases by a factor that we call Z. Namely, Qtoday = Qprm/Z. The value of Z is 
galaxy-dependent. The spherical model gives Z ~ 41000 and iV-body simulations indicate: 10000 > Z > 1 (see [1]). 
Combining the value of Qtoday and (today) with the theoretical analysis yields that m must be in the keV scale 
and Td can be larger than 100 GeV. More explicitly, we get general formulas for m and gd [1]: 

23 \/7r 1 3 _s 2t- q'i T^ i , 3 

3» gt 3s 7r2 Udm Pc 

where hn = y^" Fd{y) dy , n = 1,2 and q|„^ = zi 0.18 keV using the dSphs data, = 
0.2348 meV , ^ 0.228 and pc (2.518 meV)^. These formulas yield for rchcs decoupling UR at LTE: 

f ^ w f 0.568 3 ^1 f 155 Fermions 

" 1 5 J I 0-484 ^9d = g^Z^ 3^^^^^^ . 

Since g = 1 - 4, we see that gd > 100 ^ > 100 GeV. Moreover, 1 < < 10 for 1 < Z < 10000. For example for 
DM Majorana fermions (g = 2) m ~ 0.85 keV. 

We get results for m and gd on the same scales for DM particles decoupling UR out of thermal equilibrium [1]. For 
a specific model of sterile neutrinos where decoupling is out of thermal equilibrium: 

0.56 keV < TO^ < 1.0 kcV , 15 < < 84 

For relics decoupling non-rclativistic we obtain similar results for the DM particle mass: kcV < m < McV [1]. 
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Notice that the dark matter particle mass m and decoupling temperature Td are mildly affected by the uncertainty 
in the factor Z through a power factor 1/4 of this uncertainty, namely, by a factor 10* ~ 1.8 

The comoving free-streaming) wavelength, and the Jeans' mass are obtained in the range 

kpc < Xfsiz) < kpc , 0.45 10^ Afo < ^ < 0.45 10^ Mq . 

Vl + z + z (1 + z)+2 

These values at z = are consistent with the A^-body simulations and are of the order of the small dark matter 
structures observed today . By the beginning of the matter dominated era z ^ 3200, the masses are of the order of 
galactic masses ^ 10^^ Mq and the comoving free-streaming wavelength scale turns to be of the order of the galaxy 
sizes today ^ 100 kpc. 

Lower and upper bounds for the dark matter annihilation cross-section (Tq are derived: ctq > (0.239 — 
0.956) 10~^ GeV~^ and ctq < 3200 m GeV~^ . There is at least five orders of magnitude between them, the 
dark matter non-gravitational self-interaction is therefore negligible (consistent with structure formation and observa- 
tions, as well as by comparing X-ray, optical and lensing observations of the merging of galaxy clusters with TV-body 
simulations). 

Typical 'wimps' (weakly interacting massive particles) with mass m = 100 GeV and Td = 5 GeV would require a 
huge Z ^ 10^'^, well above the upper bounds obtained and cannot reproduce the observed galaxy properties. They 
produce an extremely short free-streaming or Jeans length Xfs today A/s(0) ^ 3.51 10~^ pc = 72.4 AU that would 
correspond to unobserved structures much smaller than the galaxy structure. Wimps result are strongly disfavoured. 

Galaxies are described by a variety of physical quantities: 

(a) Non-universal quantities: mass, size, luminosity, fraction of DM, DM core radius tq, central DM density po- 

(b) Universal quantities: surface density fiQ = rg po and DM density profiles. Mbh /Mhaio (or halo binding 
energy). The galaxy variables are related by universal empirical relations. Only one variable remains free. That is, 
the galaxies are a one parameter family of objects. The existence of such universal quantities may be explained by 
the presence of attractors in the dynamical evolution. The quantities linked to the attractor always reach the same 
value for a large variety of initial conditions. This is analogous to the universal quantities linked to fixed points in 
critical phenomena of phase transitions. The universal DM density profile in Galaxies has the scaling property: 

p{r)=poF(-) , i^(0) = l , x^-, 

where tq is the DM core radius. As empirical form of cored profiles one can take Burkert's form for F{x). Cored 
profiles do reproduce the astronomical observations. 

The surface density for dark matter (DM) halos and for luminous matter galaxies is defined as: /xqd = Poy '''o = 
halo core radius, po = central density for DM galaxies. For luminous galaxies po = Pi^o) (Donato et al. 09, Gentile 
et al. 09). Observations show an Universal value for poD- independent of the galaxy luminosity for a large number of 
galactic systems (spirals, dwarf irregular and spheroidals, elliptics) spanning over 14 magnitudes in luminosity and of 
different Hubble types. Observed values: 

POD ~ 120 ^ = 5500 (McYf = (17.6 Mev)^ , 5kpc < ro < lOOkpc . 
pc^ 

Similar values poD — 80 are observed in interstellar molecular clouds of size rg of different type and composition 

pc 

over scales 0.001 pc < rg < 100 pc (Larson laws, 1981). Notice that the surface gravity acceleration is given by poo 
times Newton's constant. 

We combine in refs.[3-4] the theoretical evolution of density fluctuations computed from flrst principles since de- 
coupling till today to the observed properties of galaxies as the surface density and core radius. We obtain that (i) 
the dark matter particle mass must be in the keV scale both for in and out of thermal equilibrium decoupling (ii) the 
density profiles are cored for keV scale DM particles and cusped for GeV scale DM particles (wimps). 
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Mass profiles and c ^/dm relation in X-ray luminous galaxy clusters 

The distribution of the total and baryonic mass in galaxy clusters is a fundamental ingredient to validate the 
scenario of structure formation in a Cold Dark Matter (CDM) Universe. Within this scenario, the massive virialized 
objects are powerful cosmological tools able to constrain the fundamental parameters of a given CDM model. The 
A'^— body simulations of structure formation in CDM models indicate that dark matter halos aggregate with a typical 
mass density profile characterized by only 2 parameters, the concentration c and the scale radius (e.g. [8], hereafter 
NFW). The product of these two quantities fixes the radius within which the mean cluster density is 200 times 
the critical value at the cluster's redshift [i.e. i?200 = C200 x Tg and the cluster's volume V — 4/37ri?2oo is equal 
to M2oo/(200/9c,2), where M200 is the cluster gravitating mass within i?20o]- With this prescription, the structural 
properties of DM halos from galaxies to galaxy clusters are dependent on the halo mass, with systems at higher 
masses less concentrated. Moreover, the concentration depends upon the assembly redshift (e.g. [1]), which happens 
to be later in cosmologies with lower matter density, firm f-^d lower normalization of the linear power spectrum on 
scale of 8h~^ Mpc, as, implying less concentrated DM halos of given mass. The concentration - mass relation, and 
its evolution in redshift, is therefore a strong prediction obtained from CDM simulations of structure formation and 
is quite sensitive to the assumed cosmological parameters (e.g. NFW; [1]; [5]). 

Recent X-ray studies (see reference in [2], [6], [9]) have shown good agreement between observational constraints at 
low redshift and theoretical expectations from Cold Darm Matter over a wide range of halo mass. For instance, it is 
confirmed at high significance that the concentration decreases with increasing mass, requiring a as, the dispersion 
of the mass fluctuation within spheres of comoving radius of 8 h^^ Mpc, in the range 0.76 — 1.07 (99% confidence) 
definitely in contrast to the lower constraints obtained, for instance, from the analysis of the WMAP 3 years data 
(see [2]). However, a possible tension between the observational constraints and the numerical predictions is also 
observed, in the sense that either the relation is steeper than previously expected or some redshift evolution has to be 
considered ([9]). Moreover, there is evidence that strong Icnsing measurements of the concentration are systematically 
larger than the ones estimated in the X-ray band, and 55 per cent higher, on average, than the rest of the cluster 
population ([3]). 

In [6] , we extend the spectral analysis presented in [7] for a sample of 44 X-ray luminous galaxy clusters located in 
the redshift range 0.1 — 0.3 with a detailed spatial analysis of the surface brightness profiles with the aim to (1) recover 
their total and gas mass profiles, (2) constraining the cosmological parameters as and fim through the analysis of the 
measured distribution of €200, M200 and baryonic mass fraction in the mass range above IO^^Mq. We note that this 
is the statistically largest sample for which this study has been carried on up to now between z = 0.1 and z = 0.3. 

We use the profiles of the spectroscopically determined ICM temperature and of the PSF-corrected surface bright- 
ness estimated to recover the X-ray gas, the dark and the total mass profiles, under the assumptions of the spherical 
geometry distribution of the Intracluster Medium (ICM) and that the hydrostatic equilibrium holds between ICM 
and the underlying gravitational potential. 

Our dataset is able to resolve the temperature profiles up to about 0.6 — 0.8i?5oo and the gas density profile, 
obtained from the geometrical deprojection of the PSF-deconvolved surface brightness, up to a median radius of 
0.9i?5oo- Beyond this radial end, our estimates arc the results of an extrapolation obtained by imposing a NFW 
profile for the total mass and different functional forms for il/gas- 

We estimate, with a relative statistical uncertainty of 15 — 25%, the concentration C200 and the mass M200 of 
the dark matter (i.e. total— gas mass) halo. We constrain the C200 ~ -^200 relation to have a normalization C15 = 

C200 X (1 + ^) X {M2oo/IO^^Mq) ^ of about 2.9 — 4.2 and a slope B between —0.3 and —0.7 (depending on the methods 
used to recover the cluster parameters and to fit the linear correlation in the logarithmic space) , with a relative error 
of about 5% and 15%, respectively. 

We put constraints on the cosmological parameters (crgi^m) by using the measurements of C200 and M200 and 
by comparing the estimated values with the predictions tuned from numerical simulations of CDM universes. In 
doing that, we propagate the statistical errors (with a relative value of about 15 — 25% at la level) and consider the 
systematic uncertainties present both in the simulated datasets (~ 20%) and in our measurements. 

When the subsample of 11 LEC clusters, that are expected to be more relaxed and with a well-formed central 
coohng core, is considered, we measure 7 = 0.56 ± 0.04, F = 0.39 ± 0.02 as = 0.83 ± 0.1 and r2„i = 0.26 ± 0.02 (at 2a 
level) . 

We show that the study of the distribution of the measurements in the c — Mdm ~ /gas plane provides a valid 
technique already mature and competitive in the present era of precision cosmology. 
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M200 (Mo) M200 (Mq) 

FIG. 5: Data in the plane (£200,^200) used to constrain the cosmological parameters (fimjCrs)- The dotted lines show the 
predicted relations from [1] for a given ACDM cosmological model at 2: = (from top to bottom: as = 0.9 and ag — 0.7). 
The shaded regions show the predictions in the redshift range 0.1 — 0.3 for an assumed cosmological model in agreement with 
WMAP-1, 5 and 3 years (from the top to the bottom, respectively) from [5]. The dashed lines indicate the best-fit range at la 
obtained for relaxed halos in a WMAP-5 years cosmology from [4; thin lines: z = 0.1, thick lines: z — 0.3]. (Left) The color 
code indicates the objects at 2 < 0.15 (blue), in the range 0.15 < 2 < 0.25 (green) and at 2 > 0.25 (red). {Right) Distribution 
of Low (LEG), Medium (MEG), High (HEC) Entropy Gore systems. Filled diamonds indicate the data where a more robust 
(i.e. with well defined and constrained free parameters) mass reconstruction is achievable. 



However, we highlight the net dependence of our results on the models adopted to relate the properties of a DM halo 
to the background cosmology. In this context, we urge the iV— body community to generate cosmological simulations 
over a large box to properly predict the expected concentration associated to the massive (> lO^^M©) DM halos 
as function of as, Q^a and redshift. The detailed analysis of the outputs of these datasets will provide the needed 
calibration to make this technique more reliable and robust. 
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The "Bosma Effect" Revisited - Correlations between the ISM and DM in Galaxies 

The problem of explaining motions within galaxies - the internal stellar dynamics of the Sun's neighborhood or 
within an elliptical galaxy and the rotation curves of spiral galaxies - was one of the original reasons for the creation 
of the cold Dark Matter paradigm, since it is on the scales of galaxies that the masses of the putative DM particles 
are most relevant. With the creation of the AC DM paradigm, much of the emphasis has shifted onto the larger 
scales of galaxy clusters and cosmological large-scale structure, but DM in galaxies is one of the critical links in the 
Borromcan chain of the Concordance Cosmological Model (CCM) of the Universe: astronomical DM; Dark Energy 
and inflationary cosmology; and the standard model of particle physics. 

Despite the astounding performance of the CCM, there are many signs that the simple picture of CDM in galaxies 
is not working. The most troubling signs of the failure of the CDM paradigm have to do with the tight coupling 
between baryonic matter and the dynamical signatures of DM in galaxies, e.g. the TuUy-Fisher relation [1], the 
stellar disc-halo conspiracy [2], the maximaum disc phenomenon [3], the MOdified Newtonian Dynamics (MOND) 
phenomenon [4], the baryonic Tully-Fisher relation [5], the baryonic mass discrepancy-acceleration relation [6], the 
1-parameter dimensionality of galaxies [7], and the presence of both a DM and a baryonic mean surface density 
[8] [9]. The strangest of these relations is the "Bosma effect" [10] [11]: the centripetal contribution of the dynamically 
insigificant interstellar medium (ISM) in spiral galaxies is directly proportional to that of the dominant DM. The 
constant of proportionality has been determined for about 100 galaxies, with dwarf galaxies showing a smaller and 
late-type spirals showing a larger factor. 

Hoekstra, van Albada & Sancisi [12] set out to test the Bosma effect, showed that it indeed allowed a very detailed 
fit to the rotation curves of many well-studied galaxies, but concluded that it was not real. Reviewing their arguments, 
it is clear that their negative judgement was very conservative - by the normal standards of rotation curves, the results 
were, in fact, very convincing. Since their fits were performed by hand and not compared with the corresponding 
CDM model fits, it was not possible to make any formal conclusions and certainly not possible to reject the effect as 
non-physical. 

Using the much better data made available by the Spitzer Infrared Nearby Galaxy Survey [13], The HI Nearby Galaxy 
Survey [14], and the analyses by de Block et al. [15], we have tested the Bosma effect and compared the results against 
standard CDM models. The use of infrared photometry and colours permits formal fits to the stellar components 
nearly independent of extinction corrections and with reasonably reliable mass-to-light estimates. In addition to the 
standard bulge, stellar disk, and visible HI components, we fitted the rotation curves with the addition of either one 
or two Bosma components, using either the HI disc (so-called "simple Bosma" models) or both the stellar and the HI 
discs ("classic Bosma" models) as proxies, where the stellar disc is used as a proxy for the molecular gas obviously 
present in regions of previous and current star formation. For comparison, we also fit the data with self-consitent 
NFW models, where the compactness of the halo is a function of the halo mass [16] or with the Burkert halo mass 
profiles used in the "Universal Rotation Curve" model [17]. 

The "simple" Bosma models using only the HI as a proxy are remarkably good in the outer discs, as shown by 
Bosma, independent of the shape of the rotation curve. However, the inner disks are not well-fit by pure "Hl-scaling" : 
the saturation of the HI surface densities above levels around 10 Mq/pc^ results in centripetal contributions which 
are clearly too small. This problem is not present in "classic" Bosma models, since the saturation of the HI profiles 
occurs exactly where the stellar component starts to dominate, permitting a perfect compensation. 

The self-consistent NFW models are clearly much worse than the Bosma models. Not unexpectedly, the fully 
unconstrained URC models perform best. The median ratio of the fitted values between the "classical" Bosma 
and URC models is only about 1.5 compared with a value of 3 between the NFW and URC fits. Thus, the "classical" 
Bosma models are only marginally worse than the best-fitting DM models. Since the Bosma model works so well 
in galaxies covering a very wide range of masses, sizes, stellar contributions, and rotation curve shapes, one must 
conclude that the Bosma effect is, in fact, real. 

At the time Bosma initially reported this effect, it was reasonable to assume that the DM might be in the discs, but 
already at the time of Hoektra, van Albada & Sancisi's attempted test, the paradigm had shifted to a quasi-spherical 
DM halo. However, the centripetal pattern of a spherical halo is dramatically different from that of a disc: in a disc, 
the matter at all radii contribute to the local centripetal effect, whereas in a spherical halo only the inner mass is felt. 
This means that a physical interpretation of the Bosma effect withing the classic DM paradigm is highly unlikely: 
how can a dominant spherical halo always teach a minor disc mass component to behave as if it were a spherical halo, 
fully independently of the properties of both. This is either a bizzare conspiracy or a potent argument against the 
DM halo paradigm. 
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If the Bosma effect is interpreted as truly being due to disc DM (presumedly baryonic, since scalable from the ISM), 
the amount of additional baryonic mass is quite significant: this implies that only 10-70% of the baryons are visible 
as stars and gas. Within the CCM, this is not a fundamental problem, since the standard paradigm has a "missing 
baryon" problem anyway [18]. The reason that a relatively modest amount of disc DM can replace a large amount 
of CDM is very simple: a disc provides much more centripetal effect per unit mass than a highly extended spherical 
halo. 

The disc DM solution has the power to explain a wide range of otherwise inexplicable phenomena, e.g. the utility 
of "maximal discs" (using the stellar disc as a proxy for disc DM is nearly paramount to fitting "maximal discs" , 
particularly since the flat HI profiles often produce a negative centripetal rotation curve contribution at small radii) , 
the extended baryonic Tully-Fisher relation, the relative value of the mean URC/DM and baryonic surface densities, 
the MOND phenomenon, the mass discrepancy-acceleration relation. Since cold gas in the ISM is very hard to detect 
- a situation which is being confirmed by new far infrared data taken by the Herschel satellite - the H2 "clumpuscule" 
model of Pfenniger, Combes & Martinet [19] provides a natural explanation for the Bosma effect. The main argument 
against disc DM is that such a massive disc must be unstable: this argument must be taken very seriously, but the 
standard arguments are not nearly as tight as is normally assumed: this will be one of the topics in a future paper 
[20]. 

The Bosma effect does not show that there is no DM, since it only probes the gravitational potentials and sources 
within galaxies. For instance, large total masses are derived for the Milky Way and M31 using the kinematics of 
satellite galaxies [21], although it is difficult to obtain tight constraints on the true total masses due to a variety of 
effects [20]. Within the Bosma effect paradigm, these large masses are most easily explained using a DM component 
whose core radii are as large or larger than the visible galaxies: DM halos with these extents have practically no 
effect on the dynamics on scales of galaxies and only show up on the scales of groups (e.g. the Local Group), clusters 
of galaxies, and cosmologically. The presentation of Moni Bidin at this conference (see [21]) fits very well with this 
model: there is no sign of DM at the radius of the Sun. 

The topic of this workshop - Warm to rather Hot DM (depending upon where one chooses to set the mass scale) - 
provides the most natural explanation for phenomena outside of the scale of visible galaxies while providing a potential 
reason for the pure dominance of baryons on the scales of visible galaxies. While the streaming scale of bare 10 keV 
neutrinos (ca. 10 kpc) is still too small, it is conceivable that processes during the formation of galaxies resulted in a 
transfer of angular momentum from the dissipative baryons to the W/HDM. 

We would like to thank the THINGS consortium, particularly E. de Blok, for making their data and rotation curve 
analyses available. 
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Light sterile neutrino as warm dark matter and the structure of galactic dark halos 

Abstract: We study the formation of nonlinear structure in A Warm Dark Matter (WDM) cosmology using large 
cosmological N-body simulations. We assume that dark matter consists of sterile neutrinos that are generated through 
nonthermal decay of singlet Higgs bosons near the Electro- Weak energy scale. Unlike conventional thermal relics, the 
nonthermal WDM has a peculiar velocity distribution, which results in a characteristic shape of the matter power 
spectrum. We perform large cosmological N-body simulations for the nonthermal WDM model. Wc compare the 
radial distribution of subhalos in a Milky Way size halo with those in a conventional thermal WDM model. The 
nonthermal WDM with mass of 1 keV predicts the radial distribution of the subhalos that is remarkably similar to 
the observed distribution of Milky Way satellites. 

Introduction 

Alternative models to the standard A Cold Dark Matter model have been suggested as a solution of the so-called 
'Small Scale Crisis'. One of them is A Warm Dark Matter cosmology, in which Dark Matter particles had non-zero 
velocity dispersions. The non-zero velocities smooth out primordial density perturbations below its free-streaming 
length of sub-galactic sizes. The formation of subgalactic structure is then suppressed. Moreover, The large phase 
space density may prevent dark matter from concentrating into galactic center. Particle physics models provide 
promising WDM candidates such as gravitinos, sterile neutrinos and so on. Gravitinos with mass of ^keV in generic 
Supergravity theory are produced in thermal bath immediately after reheating and are then decoupled kinematically 
from thermal bath similarly to the Standard Model (SM) neutrinos because gravitinos interact with SM particles only 
through gravity. The thermal relics have a Fermi-Dirac (FD) momentum distribution. Sterile neutrinos arc initially 
proposed in the see-saw mechanism to explain the masses of the SM neutrinos. If the mass of sterile neutrinos 
is in a range of ~keV and their Yukawa coupling is of order ^ 10"^*^, then it cannot be in equilibrium with SM 
particles throughout the thermal history of the universe. There are several peculiar production mechanisms for sterile 
neutrinos, such as Dodelson-Widrow (DW) mechanism [1] and EW scale Boson Decay (BD) [2]. In DW mechanism, 
sterile neutrinos arc produced through oscillations of active neutrinos, and its velocity distribution has a Fermi-Dirac 
form just like gravitinos. In the BD case, sterile neutrinos are produced via decay of singlet Higgs bosons and they 
have generally a nonthermal velocity distribution. The nonthermal velocity distribution imprints particular features 
in the transfer function of the density fluctuation power spectrum [3]. The transfer function has a cut-off at the 
corresponding free-streaming length, but it decreases somewhat slowly than thermal WDM models. In this article, 
wc study the formation of nonlinear structure for a cosmological model with nonthermal sterile neutrino WDM. 
We perform large cosmological N-body simulations. There are several models of nonthermal WDM. For example, 
gravitinos can be produced via decay of inflaton [4] or long-lived Next Lightest Supersymmetric Particle (NLSP) [5]. 
Our result can be generally applied to the formation of nonlinear structure in these models as well. 

Nonthermal sterile neutrino 

The clustering properties of the above BD sterile neutrino model is investigated by Boyanovsky [3] , who solved the 
linearized Boltzmann-Vlasov equation in the matter dominant era when WDM has already become non-relativistic. 
Firstly, solving the Boltzmann equation for sterile neutrinos produced by the singlet boson decay process, wc find the 
most of contribution to the present number of sterile neutrinos comes when the temperature decreases to the EW 
scale. We then get the velocity distribution of sterile neutrinos. The BD distribution has a distinguishable feature 
from usual Fermi-Dirac distributions at small velocities, y ~ P/T{t): 

/bd(2/) oc \, fpDiy) oc 1. 

The velocity distribution is imprinted in the power spectrum, which decreases slowly across the free-streaming length 
scale. Following Boyanovsky [3], we solve the linearlized Boltzmann-Vlasov equation. We define the comoving free- 
streaming wavenumber as 

_ r siJof^A/ "1 5 

akin to the Jeans scale at the matter-radiation equality. In the BD case, the present energy density of dark matter is 
determined by the Yukawa coupling. Assuming the relativistic degree of freedom g is the usual SM value g ~ 100 for 



17 



the m = 1 kcV sterile neutrino dark matter, we find 

fcg° = 18 /i Mpc~\ 

while for the m = 1 keV gravitino dark matter, which has the FD distribution, we should take g ~ 1000 to get the 
present energy density of dark matter. Then, for m ~ 1 keV gravitinos, 

fcf™ ^32 h Mpc~\ 

Although these two models have almost the same free-streaming length, their linear power spectra show appreciable 
differences, as seen in Fig. [51 There, we adopt the numerically fit transfer function in Bode et al. [6] for the linear 
power spectrum of the gravitino dark matter. 

The enhancement of the velocity distribution in the low velocity region leads to the slower decrease of the linear 
power spectrum. This implies that we should care not only free-streaming scale or velocity dispersion, but also the 
shape of the velocity distribution in studying the formation of the nonlinear object below the cut-off (free-streaming) 
scale. 

Using these linear power spectrum, we have performed direct numerical simulations to study observational signatures 
of the imprinted velocity distribution in the subgalactic structures. We start our simulations from a redshift of 2: = 9. 
We use N = 256'^ particles in a comoving volume of 10 Mpc on a side. The mass of a dark matter particle is 
4.53 X 10^ Mq and the gravitational softening length is 2 h~-^ kpc. 

Simulation Results 

Fig. [7] shows the projected distribution of dark matter in and round a Milky Way size halo at z = 0. The plotted 
region has a side of 2 Mpc. Dense regions appear bright. Clearly, there are less subgalactic structures for the 
WDM models. The 'colder' property of the nonthermal WDM shown in the linear power spectrum (see Fig. ^ can 
also be seen in the abundance of subhalos. 

In Fig. m we compare the cumulative radial distribution of the subhalos in our 'Milky Way' halo at 2; = with 
the distribution of the observed Milky Way satellites [7]. Interestingly, the nonthermal WDM model reproduces 
the radial distribution of the observed Milky Way satellites in the range above ^ 40 kpc. Contrastingly, the CDM 
model overpredicts the number of subhalos by a factor of 5 than the observed Milky Way satellites. This is another 
manifestation of the so-called 'Missing satellites problem'. The thermal WDM model surpresses subgalactic structures 
perhaps too much, by a factor of 2 — 4 than the observation. 

Summary 

We have studied the formation of the nonlinear structures in a AWDM cosmology using large cosmological N-body 
simulations. We adopt the sterile neutrino dark matter produced via the decay of singlet Higgs bosons with a mass 
of EW scale. The sterile neutrinos have a nonthermal velocity distribution, unlike the usual Fermi-Dirac distribution. 
The distribution is a little skewed to low velocities. The corresponding linear matter power spectrum decreases slowly 
across the cut-off scale compared to the thermal WDM, such as gravitino dark matter. Both of the two models have 
the same mass of 1 keV and an approximately the same cut-off (free-streaming) scale of a few 10 h Mpc~^. We 
have shown that this 'colder' property of the nonthermal WDM can be seen in richer subgalactic structures. The 
nonthermal WDM model with mass of 1 keV appears to reproduce the radial distribution of the observed Milky Way 
satellites. 
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FIG. 6: The linear power spectra for the nonthermal WDM (red line), thermal WDM (green line) and CDM (blue line). The 
dark matter mass of the both WDM model is m = 1 keV. 




FIG. 7: The projected distributions of the substructures in our 'Milky Way' halo at 2; = 0. The sidelength of the shown region 
is 2 Mpc. The results are for three dark matter models, CDM, nonthermal WDM with 1 keV mass and thermal WDM 
with IkeV mass, respectively, from left to right. 
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FIG. 8: The radial distribution of the subhalos in our 'Milky Way' halo at z = Q. For simulation results for CDM (black), 
IkeV nonthermal WDM (red), and IkeV thermal WDM (blue). Green bars show the distribution of the observed Milky Way 
satellites. 
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Larson' laws and the universality of molecular cloud structures 

It has long been recognized that star formation is inextricably linked to the molecular clouds where the process is 
taking place, and therefore it is important to study the structure of these objects. One of the first attempts in this 
direction has been carried out by Larson 1981. In his seminal work, Larson used molecular line data available from 
earlier studies (mostly millimeter data of nearly objects) and showed that molecular clouds obey three scaling relations: 
(1) a power-law relationship between the length L of the cloud and its velocity dispersion CTv, with CTv oc L^-^^; (2) 
approximate virial equilibrium, with 2GM j u^L ~ 1; and (3) a relationship between the density n of the cloud and its 
length, with n oc i^^-^. Larson's third law, which is the main focus of this talk, implies that molecular clouds have 
approximately constant column densities E, since E ~ nL cx L~^'^ . 

Since their formulation, Larson's laws have been the subject of several observational and theoretical studies. From 
the observational point of view, Solomon et al. 1987 presented ^^CO data for a 273 nearby clouds from the University 
of Massachusetts-Stony Brook (UMSB) Galactic Plane Survey Sanders et al. 1986. They found a size-line width 
relationship with a power index (0.5) steeper than the one derived by Larson 1981. Additionally, in agreement with 
Larson's third law, they found that the molecular gas surface density is approximately constant for all clouds with 
S(H2) = 170 M0 pc~^. Recently, the same sample of clouds has been reanalysed by Heyer et al. 2009 using data 
from the Boston University-FCRAO Galactic Ring Survey Jackson et al. 2006. The use of ^^CO (J — 10) emission 
instead of ^^CO ensures that a large fraction of the data are optically thin; additionally, the data used have a much 
higher angular sampling and spectral resolution. Heyer et al. 2009 confirmed Larson's third law with a relative scatter 
(approximately a factor 3) similar to previous studies. However, surprisingly they found a median mass surface density 
of molecular hydrogen for this sample of 42 M© pc^, thus significantly smaller than the one derived by Solomon et 
al. 1987. 

On the theoretical side, there have been many attempts to explain Larson's laws using numerical simulations. In 
many cases, the validity of Larsons relations, and especially of the third law, has been questioned Kegel et al. 1989 
- Ballesteros-Paredes et al.2006. In particular, it has been suggested that this law is merely the result of the limited 
dynamic range of observations, and that in reality mass surface densities of molecular clouds span at least two orders 
of magnitude. 

We re-examine the validity of Larson's third law using extinction as a tracer of molecular gas Lada et al. 1994. 
The use of this tracer, in combination with advanced techniques Lombardi et al. 2001 and 2009, allows us to probe 
clouds over a large dynamical range (typically more than two order of magnitudes in extinction); additionally, the 
column density measurements use a simple and reliable tracer, dust. 

We consider first the following version of Larson's third law. Since we have at our disposal complete extinction 
maps, we can consider the area 5 of a cloud above a given extinction threshold Aq (unless otherwise noted, we will 
refer to extinction measurements in the K band, Ak, and drop everywhere the index K). We then define the cloud 
size implicitly from S = 7r(i/2)^ (or the cloud radius as i? = L/2). Similarly, we can consider the cloud mass M 
above the same extinction threshold. 

Figure [S] left shows the amount of mass different clouds have above extinction thresholds of Ak = 0.1 mag and 
Ak = 0.5 mag as a function of the cloud "radii", together with the best power-law fit. As apparent from this plot, 
all clouds follow exquisitely well a Larson- type relationship, with Al (x R^, and have therefore very similar projected 
mass densities at each extinction threshold. The exceptionally small scatter observed in Fig. [S] is also confirmed by a 
quantitative analysis: at all extinctions considered, data follow the best-fit power-laws with relative standard deviations 
always below 15%. 

Figure ini right shows the second version of Larson's third law considered here, i.e. the mass vs. radius relationship. 
As apparent from this figure, the tracks for the various clouds have similar trends, but span a relatively large range 
of masses. In the range R E [0.1, 1] pc we can fit a power-law of the form M{R) = 380 Mq (i?/pc)^'^, a result that 
compares well with the one obtained by Kauffmann et al. 2010, M{R) = 400 M© (-R/pc)^-^. Different clouds have 
quite similar exponents (the standard deviation of the power-law index is 0.18), but rather different masses (the 
best- fit scale parameter for the mass ranges from 150 to 710 M©). Note, however, that since the power-law index is 
significantly different from two, errors on the assumed distances of the clouds would affect the scale parameter for the 
mass. From this analysis we conclude that Larson's third law is not an accurate description of the mass vs. radius 
relationship for single clouds. Specifically, at larger scales all clouds show a fiattening of the curves and deviates 
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significantly from a power-law, while at smaller scales clouds follow power-laws, but with an exponent significantly 
different than two. 

Our results can be summarized as follow: 

1. Using near- infrared extinction maps of a set of nearby clouds we tested Larson's third law for molecular clouds, 
the constancy of average mass surface densities above a given extinction threshold. We verified this scaling 
law to a relatively high degree of precision. We found a very small (< 15%) relative scatter for the measured 
column densities independent of the adopted extinction thresholds over a very large range, from Ak = 0.1 mag 
to Ak = 1.5 mag. Additionally, we found the value of the average mass surface density to be a function of the 
adopted extinction threshold. 

2. We verified that Larson's third law does not hold when considering the mass-radius relation within single clouds. 
In the range R G [0.1, 1] pc we find that the mass scales as M{R) cx R^^^ , and is therefore significantly shallower 
than what was predicted by Larson; at larger radii, the relation appears to fiatten even more. 

3. We interpret these results, and in particular item 1 above, as the effects of a universal physical structure 
shared among the different clouds. This universal structure is represented by a uniformity in the cloud density 
distributions. We find that a log-normal model is able to account for this uniformity, provided that the log- 
normal parameters are restricted to relatively narrow ranges. This suggests that Larson's third law might be a 
consequence of this special property of cloud structure. 
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Warm Dark Matter with future cosmic shear data 

Free-streaming dark matter particles dampen the overdensities on small scales of the initial linear matter density 
field. This corresponds to a suppression of power in the linear matter power spectrum and can be modelled relatively 
straightforwardly for an early decoupled thermal relic dark matter particle. Such a particle would be neutrino-like, 
but heavier; an example being the gravitino in the scenario, where it is the Lightest Supersymmetric Particle and 
it decouples much before neutrinos, but while still relativistic, see [1]. Such a particle is not classified as Hot Dark 
Matter, like neutrinos, because it only affects small scales as opposed to causing a suppression at all scales. However 
its free-streaming prevents the smallest structures from gravitationally collapsing and does therefore not correspond to 
Cold Dark Matter. The effect of this Warm Dark Matter may be observable in the statistical properties of cosmological 
Large Scale Structure. 

Cosmic shear is the weak gravitational lensing of the images of very distant galaxies caused by the deflection of 
photons from these galaxies by the gravitational potential wells of the intervening dark matter density field. Because 
this effect does not strongly depend on baryonic physics it is a promising probe of the statistics of the dark matter 
density field. Future cosmic shear data will be able to measure the cosmic shear power spectrum to a very good 
accuracy and will therefore provide a probe for the non-linear matter power spectrum at low redshifts (late times). 
Unfortunately robust modelling of gravitational non-linearities at late times is not straightforward, but is absolutely 
necessary in order to extract cosmological information from large scale structure probes, in particular from cosmic 
shear. 

The suppression of the linear matter density field at high redshifts in the WDM scenario can be calculated by solving 
the Boltzmann equations. A fit to the resulting linear matter power spectrum in the simple thermal relic scenario is 
provided by [2]. This linear matter power spectrum must then be corrected for late-time non-linear collapse. This 
can be done using the halofit method of [3] or using the halo model ([4]). In [5], both of these non-linear methods are 
used to calculate weak lensing effect from a future survey and a limit of towdm ^ 2 keV is predicted. WDM with 
greater particle masses in the simple thermal relic scenario is predicted not to be detectable with the next generation 
of weak lensing surveys. 

However, both of the above non-linear methods were developed assuming CDM and are therefore not necessarily 
appropriate for the WDM case. For this reason, [6] modify the halo model. Firstly, they treat the dark matter density 
field as made up of two components: a smooth, linear component and a non-linear component, both with power at 
all scales. Secondly, they introduce a cut-off mass scale, below which no haloes are found. Thirdly, they suppress the 
mass function also above the cut-off scale and finally, they suppress the centres of halo density profiles by convolving 
them with a Gaussian function, whose width depends on the WDM relic thermal velocity. The latter effect is shown 
to not be significant in the WDM scenario for the calculation of the non-linear matter power spectrum at the scales 
relevant to weak lensing. They make predictions for the next generation weak lensing surveys, which are consistent 
with the results of [5]. See Figure [TU] for the plot of the weak lensing power spectra in different WDM scenarios. 

In order to determine the validity of the different non-linear WDM models, [7] ran cosmological simulations with 
WDM. They provide a fitting function that can be easily applied to approximate the non-linear WDM power spectrum 
at redshifts z = 0.5 — 3.0 at a range of scales relevant to the weak lensing power spectrum. Figure [TT] shows the 
percentage differences between the WDM and CDM non-linear matter power spectrum in the simple thermal relic 
scenario for different WDM masses. 
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FIG. 10: Comparison of the weak lensing convergence power spectrum in the CDM and WDM halo models, as a function of 
angular multipole for source galaxies at = 1. Top panel: Absolute power. Dashed and solid lines show the predictions from 
linear theory and from the nonlinear halo model. Black lines represent results for CDM. Coloured lines with red, green, blue, 
cyan, magenta denote results for the WDM model with particle masses mwDM = 0.25,0.5,0.75,1.0,1.25 keV respectively. 
Bottom panel: Ratio of the WDM model predictions to the CDM predictions. Lines styles and colours are as above. 
This figure is from [6]. 
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FIG. 11: Percentage difference between warm dark matter non-linear power and cold dark matter for the different runs. The 
resolution is kept fixed in this plot and only 25 Mpc/h boxes are considered. Orange, green, black, blue and red curves refer 
to mwDM = 0.25, 0.5, 1, 2, 4 keV, respectively. The continuous lines represent the large scale estimate of the power, while the 
dashed ones describe the small scale power obtained with the folding method. The four panels represent different redshifts at 
2 = 0, 1, 2, 5 (bottom right, bottom left, top right and top left, respectively). The dotted coloured curves plotted at 2 = and 
2 = 5 are the redshift independent linear suppression between the different models. 
This figure is from [7]. 
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Neutrino Model Building and keV sterile Neutrino Dark Matter 

One of the most exciting problems in astrophysics and cosmology, but maybe also one of the most attractive 
connection points of these fields to elementary particle physics is the explanation of the so-called Dark Matter (DM) 
in the Universe. DM could be cold (non-relativistic), like in the standard cosmological model, or hot (relativistic) , 
which is, however, strongly disfavored by structure formation. Between those limiting cases there would be warm 
Dark Matter (WDM), favored by simulations of large scale structure formation. A particularly interesting candidate 
for WDM from the particle physics side is a sterile neutrino with a mass of a few keV, which can be accommodated 
for by a simple extension of the Standard Model (SM) of Elementary Particle Physics dubbed vMSM [1]. 

Indeed, gauge extensions of the SM could easily include such a scenario [2], while still being consistent with all 
phenomenological, observational, and experimental constraints. Although these frameworks can be brought into 
accordance with all data, they actually do not give an "explanation" for the values of the parameters involved, and 
actually they are scenarios rather than models: In a scenario, all the parameters involved can assume values such 
that there is no disagreement with any observation or experiment, while models give, at least to some extend, reasons 
for why certain parameters have their respective values. On the other hand, the experimental data from neutrino 
experiments seems to indicate that there is a reason for the parameter values observed [3], which makes it tempting 
to develop dedicated models with certain features that enable them to relate the different observables in such a way 
as to explain part or all of the data. 

Although the neutrino model building industry is very active [4], and although sterile (right-handed) neutrinos 
most probably exist anyway due to the necessity of a non-zero active neutrino mass, up to now, there are only a few 
models around that can explain such a peculiar mass patterns as needed for keV sterile neutrino DM. Observational 
constraints require one neutrino to be at the keV-scale, while the other two sterile neutrinos (in a three-generation 
model) must have masses of at least GeV [2]. From a theory point of view, one would expect them to be even 
heavier [4]. 

The first class of models yielding an explanation of this mass pattern is the one proposed by Kusenko et al. [5] , who 
used the exponential warping factor present in Randall-Sundrum models of extra dimensions to suppress the lightest 
sterile neutrino mass down to the keV scale, while the other two can have masses of about 10^^ GeV or higher. This 
leads to a generation dependent exponent a^, which suppresses the natural mass scale Mg according to 



such that a mild hierarchy ai > a2 > of the exponents will translate into a very strong hierarchy Mi ^ M2 ^ M3 
of the 4D sterile neutrino masses. A further bonus of this model is that, by the structure of the exponential suppressions 
involved, the seesaw mechanism to explain the light neutrino masses is not spoiled by the presence of the keV-scale 
sterile neutrino. 

A second class yielding an explanation for the pattern required has been constructed by Lindner et al. [6], who 
used an = Le — Lfj^ — Lr flavour symmetry as basic ingredient. The mechanism is illustrated in Fig. 112b .: The 
exact J-' symmetry (black) leads to a heavy (and light) neutrino mass pattern that is already quite close to the one 
required, by having one neutrino exactly massless while the other two arc degenerate and have, let us say, the mass M. 
However, instead of one sterile neutrino being massless, one would require it to be at the keV scale. Furthermore, light 
neutrino data excludes the existence of two strictly degenerate light neutrinos. Accordingly, as practically any known 
flavour symmetry, the J- symmetry has to be broken at some level. A convenient way to parametrize such a breaking 
is by introducing so-called soft breaking terms, whose magnitude s is much smaller than the magnitude M ! of the 
symmetry-preserving terms. One could also try to construct an explicit scalar sector that leads to the correct breaking, 
but the resulting neutrino mass matrix will just be an explicit realization of the soft-breaking parametrization. They 
key point is that the terms of order s will slightly alter the mass eigenvalues, in such a way that the degeneracies are 
broken and the otherwise massless neutrino obtains a small mass of magnitude s <^ M . The parameter s can easily 
be taken to have a value of a few keV, while M must be considerably larger, giving a motivation to have it at the 
GeV-scale or higher. 

This situation is depicted by the red part of Fig. [T^ . which shows the modification introduced by the broken 
symmetry. Indeed, the sterile neutrino mass eigenvalues are changed according to 



(0, M, M) ^ (s, M - s, M + s). 
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FIG. 12: The mass shifting schemes for the Le ~ — Lt and the Froggatt-Nielsen models. 



with an analogous modification in the hght neutrino sector. This model explains two of the three neutrino mixing 
angles in a natural way, by predicting ^13 = and ^23 = 45°, while unfortunately also predicting 9i2 = 45°, which is 
off compared to experiments by 6cr. This problem, however, can be cured by slightly modifying the charged lepton 
matrices, which are naturally expected to receive corrections from renormalization group running. Summing up, this 
model cannot only explain the sterile neutrino mass pattern, but it can also naturally explain part of the Icptonic 
mixing, and it docs in fact predict an exact light neutrino mass pattern, to be probed within the next few years. 

Finally, a third class of models has just been developed [7], which is based on the Froggatt-Nielsen (FN) mech- 
anism [8]. The model makes use of the generation dependent FN-suppressions to push one sterile neutrino mass 
down to the keV scale, as displayed in Fig. [I2b . The key point is that a suitable choice of FN-charges leads to a 
suppression of certain mass matrix elements by roughly one order of magnitude per charge unit. In order to achieve 
the required sterile neutrino splitting, it is necessary for one eigenvalue of the right-handed neutrino mass matrix to 
have an FN-charge that is larger than the one of the second-to lightest eigenvalue by at least six units. 

Although models based on the FN-mechanism often suffer from a high degree of arbitrariness, and hence from 
a lack of predictivity, one can show that, for the keV sterile neutrino DM framework, the number of possibilities 
can be reduced considerably by carefully taking into account all constraints. For example, in the framework under 
discussion, the FN-mechanism is incompatible with Left- Right symmetry, which has been used in the scenario of [2], 
it is in conflict with the bimaximal mixing predicted by the model from Ref. [6] , and it favors Grand Unified Theories 
based on SU (5) compared to the ones based on 50(10). One further feature is that, due to the construction principle 
of the FN-mechanism, also in this model the validity of the seesaw mechanism is guaranteed, even though one could 
naively expect that a keV sterile neutrino could spoil its success. 

In conclusion. kcV sterile neutrinos as WDM allow for very interesting connections between the Dark Matter 
problem and neutrino model building. In particular, the strong entanglement of the light neutrino data with the 
properties of the sterile neutrinos provides an interesting possibility to test models with keV sterile neutrino DM 
without the need to rely on direct (and rare) signals from the sterile neutrinos themselves. Currently, there are only 
three working models around that are able to explain the distinctive features of the keV sterile neutrino DM scenario, 
and the community should hunt for more such models in the future. 
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No evidence of dark matter in the Galactic disk 

Measuring the matter density of the Galactic disk in the proximity of the Sun through stellar kinematics is an old 
art, dating nearly a century [1,2]. The comparison of the results with the expected amount of visible matter can 
provide an estimate of the dark matter (DM) density in the analyzed volume. So far, all but two estimates [3,4] have 
converged to the conclusion that "there is no evidence for a significant amount of disk DM" [e.g., 5,6]. Despite the 
general agreement, this result is still very general, and only a couple of works obtained some stronger constraints on the 
fundamental properties of the DM halo, such as its flattening and local density [7,8]. However, previous investigations 
did not directly measure the mass density, but this was derived adjusting the Galactic model parameters to fit the 
observations, with the exception of [9]. Moreover, the results are mainly limited by the approximations used in the 
calculations, whose validity breaks down at increasing height from the Galactic plane (z), and no study has ever 
exceeded z=l.l kpc. In this restricted volume the expected amount of DM is small, and within errors the presence of 
a classical DM halo is often not required, but always allowed. 

We propose a new formulation for the estimate of the disk mass density, based on the distribution of stars in the six- 
dimentional phase space. By means of the Poisson and Jeans equation, and very basic assumptions both theoretically 
and observationally confirmed, we obtain the exact (not approximated) analytical expression of the surface mass 
density S (M© pc^^) valid at any z. The information required to solve this equation is the vertical trend of the three 
components of the dispersion ellipsoid in cylindrical coordinates ajj^z), (Ty(z), aiyiz), and of the non-diagonal term 
(T^^(z), the radial derivative of ajj, cry, and (jfj^, plus three parameters: the solar Galactocentric distance Rq. and 
the scale length and height of the disk population under study (h/{_p and h^.p, respectively). 

We appHed our formulation to a sample of thick disk red giant stars with 2MASS photometry [10], SPM3 absolute 
proper motion [11], and radial velocity measurements [12]. We thus measured S(z), for the first time, up to 4 kpc 
from the plane. The data provided no information about the radial behavior of the dispersions, but the assumed 
radial decay is supported by both observations [13] and theoretical models [14]. The results were first presented in 
[15], and a publication with all the details of the analysis is in preparation. 

As shown in Figure [T51 we find a striking coincidence between the calculated curve and the expectations for visible 
mass alone. A classical DM halo is excluded at a Qa level, and the amount of local DM allowed by the la error is 
negligible (<1 Mq pc~"^). We find that the results can be forced to allow a certain amount of local DM. but only under 
a series non-standard assumptions, e.g. a very thin and extended stellar thick disk (h^j p >4.5 kpc, h^ p <0.6 kpc) 
ruled out by observations. While the break-down of a single assumption is possible, the required series of ad-hoc 
alternative hypothesis is very unlikely. 

Our work, while representing a noticeable improvement to the knowledge of the Galactic mass spatial distribution, 
still lacks the information required to solve the equations with observed quantities only. Present and future extensive 
surveys such as SDSS [16], RAVE [17], and GAIA, will provide the missing data, allowing the study of the density 
distribution of the DM halo with unprecedented detail. 



27 



150 



100 - 



o 



N 



50 - 




VISIBLE MASS 



_] I L L 



J I I L. 



_] I I L. 



2 3 
z [kpc] 



FIG. 13: Measured surface mass density as a function of distance from the plane (full black curve). The la and 3a strips are 
indicated by the dashed and dotted curves, respectively. The lower grey curve shows the expectations for visible mass alone, 
while the upper one corresponds to a model where a classical DM halo from [18] is added. More details about these models 
can be found in [15]. 
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The MARE experiment and its capabilities to measure the mass of light (active) and heavy (sterile) 

neutrinos 

Direct neutrino mass measurement analyzing the kinematics of electrons emitted in beta decays is the most sensitive 
model independent method to assess the neutrino mass absolute value. In practice this method consists in measuring 
the minimum energy carried away by the anti-neutrino - i.e. its rest mass - by observing the highest energy electrons 
emitted in the decay. 

To date, the most sensitive experiments were carried out analyzing the decay in magnetic adiabatic collimation 
spectrometers with electrostatic filter, yielding an upper limit on the electron anti-ncutrino mass of 2.2 eV. Starting 
from 2013 or 2014 the new experiment KATRIN will analyze the "^H beta decay end-point with a much more sensitive 
electrostatic spectrometer and with an expected statistical sensitivity of about 0.2 eV. However a severe limitation to 
these experiments resides in their configuration whereas the "^H source is external to the spectrometer. This causes 
many systematic uncertainties because the measured electron energy has to be corrected for the energy lost in exciting 
atomic and molecular states, in crossing the source, in scattering through the spectrometer, and more. Therefore, to 
improve the sensitivity, it is necessary to reduce both the systematic and the statistical uncertainties. Because of the 
large weight of systcmatics, it is inherent in this type of measurement that confidence in the results can be obtained 
only through confirmation by independent experiments. 

An alternative approach is calorimetry where the beta source is embedded in the detector so that all the energy 
emitted in the decay is measured, except for that taken away by the neutrino. In this configuration, the systematic 
uncertainties arising from the electron source being external to the detector are eliminated. On the other hand, since 
calorimeters detects all the decays occurring over the entire beta energy spectrum, the source activity must be limited 
to avoid spectral distortions and background at the end-point due to pulse pile-up. As a consequence the statistics 
near the end-point is limited as well. Since the fraction of decays in a given energy interval /S.E below the end-point 
Q is only oc {/S.E/QY, the limitation on the statistics may be partially overcome by using as beta source ^^^Re, the 
beta-active nuclide with the second lowest known transition energy [Q ^ 2.5 keV). 

A perfect practical way to make a calorimetric measurement is to use thermal detectors. At thermal equilibrium, 
the temperature rise of the detector - measured by a suitable thermometer ~ is due to the sum of the energy of the 
emitted electron and of all other initial excitations. The measurement is then free from the systematics induced by 
any possible energy loss in the source and is not affected by problems related to decays on excited final states. 

Thermal detectors with absorbers containing natural Rhenium are the most straightforward way to make a calori- 
metric measurement. ^^^Rc natural isotopic abundance (A.I.= 63%) and its half-life time (ti/2 = 43.2 Gy) make it 
in principle perfect to design small sized (~ mg) high performance detectors with an activity Ap of few decays per 
second. The high statistics required can be accumulated with the help of large arrays of thermal detectors. 

To date, only two experiments have been carried out with thermal detectors containing ^*''Re: the MANU [1] and 
MIBETA [2] experiments. The two experiments, with a statistics corresponding to about lO*" decays, yielded limits 
on of about 26 cV at 95% CL and 15 eV at 90% CL respectively. The systematics affecting these experiments are 
still small compared with the statistical errors. The main sources are the background, the detector response function, 
the theoretical spectral shape of the ^^^Rc /3 decay, the Beta Environmental Fine Structure (BEFS), and the pile- up. 

The Microcalorimeter Arrays for a Neutrino Mass Experiment (MARE) project was launched by a large international 
collaboration in 2005 with the aim of measuring the neutrino mass with a calorimetric approach [3] . The baseline of 
the MARE project consists in a large array of Rhenium based thermal detectors, but different options for the isotope 
are also being considered [5]. Left panel of Fig. 1141 shows what are the experimental requirements for this experiment 
in terms of total statistics N^v, energy resolution AE and pile-up level fpp ~ i.e. the fraction of pile-up events with 
respect to the total decays given by trA^, where tr is approximately the signal rise time [4]. It is apparent that, in 
order to reach sub-eV sensitivities, the MARE experiment must collect more than 10^^ events. This calls for the use 
of large arrays - for a total number of channels of the order of 10000 - coupled to an appropriate signal multiplexing 
scheme to avoid running into insurmountable cryogenic and economic problems. Altogether these specifications are 
about the same as for IXO - the next generation X-ray space observatory - and largely proved to be technically 
feasible. The major open issue remains the metallic Rhenium absorber coupling to the sensor. In fact, in spite 
of the many efforts, persisting technical difficulties prevented from realizing the target performances with Rhenium 
absorbers so far. In order to have a viable alternative to the baseline MARE design using Rhenium (3 decay, the 
MARE collaboration is considering the possibility to use of ^^"^Ho electron capture (EC). ^^•^Ho decays to ^^'^Dy with 
a half life of about 4570 years and, because of the low transition energy, capture is only allowed from the M shell or 
higher. The EC decay may be detected only through the mostly non-radiative atomic de-excitation of the Dy atom 
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FIG. 14: Left: statistical sensitivity of a ^*^Re neutrino mass calorimetric experiment estimated with both Montecarlo (points) 
and analytic (curves) approaches. Right: statistical sensitivity to the emission of heavy neutrinos with mass mn. From lower to 
upper curve: Montecarlo symulation with Ne^ — 10^'', fpp — 10^^ and A_E = 1.5 eV; Montecarlo symulation with A'^^^ = 8x 10^, 
fpp = lO"'^ and A.E = 15 eV; MIBETA experiment; MANU experiment [2]. 



and from the Inner Bremsstrahlung (IB) radiation. It has been proposed in [6] that, thanks to its very low transition 
energy {Qec ~ 2.5 keV), the analysis of the end-point of the calorimetric spectrum can yield a very high sensitivity 
to the neutrino mass. As for /3 decay end-point experiments, also in ^^'^Ho experiments the sensitivity on the neutrino 
mass depends on the fraction of events at the end-point which increases for decreasing Qec- Unfortunately, Qec is 
only approximately known: its determinations span from 2.2 to 2.8 keV with a recommended value of about 2.555 kcV. 
In case of the more favorable values of Qec, ^^'^Ho represents a very interesting alternative to ^^^Rc. Because of the 
relatively short half life, detectors may be realized by introducing only few ^^^Ho nuclei (about 10^^ for 1 decay/s) in 
low temperature microcalorimeters optimized for low energy X-ray spectroscopy, without any further modification. 

The calorimetric spectra measured by MARE with either ^^''Re or ^^'^Ho are also suitable to investigate the emission 
of heavy (sterile) neutrinos with a mixing angle 9. Assuming the electron neutrino i^e is a mixture of two mass 
eigenstatcs i^h and i^^. with masses rriH S> m^, then ly^, ~ vj^cosO + i/HsinO and the measured energy spectrum is 
N{E, niL, niH, 6) = N{E, m^) cos^ 9 + N{E, mu) siv? 9. The emission of heavy neutrinos would manifest as a kink in 
the spectrum at an energy of Q — mn for heavy neutrinos with masses between about and Q — Eth <~ 2.5 kcV, 
where Eth is the experimental energy threshold. Right panel in Fig. [M] shows the statistical sensibility achievable by 
the MARE experiment with ^®^Re as estimated by Montecarlo simulations. Similar curves can be obtained for ^^■^Ho. 
The two curves in the upper right corner of Fig. [H are the experimental limits set by the MANU [2] and MIBETA 
experiments. The search for heavy neutrinos may be affected by systematic uncertainties due to the background and 
to the ripple observed in the ^^''Re spectrum and caused by the BEFS. 

The MARE project is subdivided into two phases. The second phase - MARE-2 - is the final large scale experiment 
with sub-electron sensitivity. The first one - MARE-1 - is a collection of activities with the aim of sorting out both 
the best isotope and the most suited detector technology to be used for the final experiment [5] . On the short term, 
intermediate scale experiments are being prepared, with the purpose of both reaching a neutrino mass sensitivity of 
few electronvolts and studying in depth the possible sources of systematic uncertainties. The first of these experiments 
is starting in Milano using two 6x6 silicon implanted arrays with AgRe04 absorbers. The experiment may be later 
expanded to up to 8 arrays, for a total of 288 crystals. The potential statistical sensitivity to heavy neutrinos for the 
full scale experiment is shown by the cross-hatched curve in Fig. [U (right panel) . 
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The effects of free streaming on warm dark matter haloes: a test of the Gunn-Tremaine limit 

The free streaming of warm dark matter particles dampens the fluctuation spectrum, flattening the mass function of 
haloes and imprinting a fine grained phase density (PSD) limit for dark matter structures. The Gunn-Tremaine limit 
is expected to imprint a constant density core at the halo center. In a purely cold dark matter model the fine grained 
phase space density is effectively infinite in the initial conditions and would therefore be infinite everywhere today. 
However when the phase space density profiles are computed using coarse grained averages that can be measured from 
N-body simulations, the value is finite everywhere and even falls with radius with a universal power law slope within 
virialised structures (Taylor & Navarro 2001). The coarse grained phase space is an average over mixed regions of 
fine grained phase space, so this behaviour is as expected (Tremaine & Gunn 1979). 

Using high resolution simulations of structure formation in a warm dark matter universe (movies available: 
|http: / /ob swww.unige.ch/^paduroiui ) we explore these effects on structure formation and the properties of warm 
dark matter halos. 

• The finite initial fine grained PSD is a also a maximum of the coarse grained PSD. resulting in PSD profiles of 
WDM haloes that are similar to CDM haloes in the outer regions, however they turn over to a constant value 
set by the initial conditions 

• The turn over in PSD results in a constant density core with characteristic size that is in agreement with the 
simplest expectations 

• We demonstrated that if the primordial velocities are large enough to produce a significant core in dwarf galaxies 
i.e. ^ kpc, then the free streaming erases all perturbations on that scale and the haloes cannot form 

• Halo formation occurs top down on all scales with the most massive haloes collapsing first 

• The concentration - mass relation for WDM haloes is reversed with respect to that found for CDM 

• Warm dark matter haloes contain visible caustics and shells 

We ran two suites of simulations, first with a IGO'^ particles in 40 Mpc box, and the second one with 300"^ in a 
42.51 Mpc box. We adopt a flat ACDM cosmology with parameters from the first year WMAP results (Spcrgcl et al. 
2003). The transfer function in WDM model has been computed using the fitting formula suggested by Bode, Turok 
and Ostriker (2001) where a, the scale of the break, is a function of the WDM parameters (Viel et al (2005)), while 
the index v is fixed. From these simulations several galaxy mass haloes were re-simulated at higher resolution, with 
and without thermal velocities. The density profile for a 7 x IO^'^^Mq is shown in Figure [15] in four different contexts, 
the CDM case, the WDM case with just the cutoff in the power spectrum as expected for a 200eV particle (WDMl), 
with velocities corresponding to the 200eV particle (WDM2), and with velocities artificially increased such that they 
correspond to a 20eV particle but with the power spectrum of the 200eV case (WDM3). 

Figure [TCI shows the corresponding coarse grained PSD profiles calculated by spherical averaging the quantity p/a^ . 

In the case of our simulations, for a constant density in the initial conditions, the phase space density is: 
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For a critical density pcrit 



1.4 X 10 MqIpc" and i7 = 0.268 we find the phase space density: 
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FIG. 15: The spherically averaged density profiles for CDM, WDMl, WDM2 & WDM3 haloes. The resolution limit is at 
approximately 0.5% of the virial radius (the softening radii are a 0.26% of the virialized radius). 
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FIG. 16: "Phase-space density" (PSD) profiles of the same haloes shown in Figure [15] calculated using p/a^. 

Thus for a particle with rrix ~ 20 cV and a ~ 150 km/s the phase space density is 2.4 x lO~^M0pc~'^ (km/s)~^ 
which gives a theoretical value for the core radius of 9.3 kpc. If we consider the core radius to be the radius where 
the density starts decreasing from the constant value, a "by eye" fit gives a slightly larger radius of 6.5 kpc, while 
the radius where the density drops by a factor of 2 is 9.4 kpc in agreement with the theoretical predictions. 

References 

1 Bode, P., Ostriker, J. P., & Turok, N. 2001, Astrophys. J. , 556, 93 

2 Dalcanton, J. J. & Hogan, C. J. 2001, Astrophys. J. , 561, 35 

3 Navarro, J. F., Frenk, C. S., & White, S. D. M. 1996, Astrophys. J. , 462, 563 

4 Strigari, L. E., Bullock, J. S., Kaplinghat, M., Kravtsov, A. V., Gnedin, O. Y., Abazajian, K. & Klypin, A. A. 

2006, Astrophys. J. , 652, 306 

5 Taylor, J. E. & Navarro, J. F. 2001, Astrophys. J. , 563, 483 

6 Tremaine, S. & Gunn, J. E. 1979, Phys.Rev.Let., 42, 407 

7 Viel, M., Lesgourgues, J., Haehnelt, M. G., Matarrese, S., & Riotto, A. 2006, Physical Review Letters, 97, 071301 

8 Wang, J., White, S. D. M. 2007 (astroph-0702575) 



32 



M. Henri Plana 

Laboratorio de Astrofi'sica Observacional Teorica e Observacional 
Univcrsidadc Estadual de Santa Cruz - Ilheus - Brazil 

Mass distribution of GalsLxies in Hickson Compact Groups 

This study presents the mass distribution for a sample of 18 late-type galaxies in nine Hickson Compact Groups [1]. 
We used rotation curves from high resolution 2D velocity fields of Fabry-Perot observations and J-band photometry 
from the 2MASS survey, in order to determine the dark halo and the visible matter distributions. The study compares 
two halo density profile, an isothermal core-like distribution and a cuspy one. We also compare their visible and dark 
matter distributions with those of galaxies belonging to cluster and field galaxies coming from two samples: 40 cluster 
galaxies of [1] and 35 field galaxies of [3]. The central halo surface density is found to be constant with respect to 
the total absolute magnitude similar to what is found for the isolated galaxies. This suggests that the halo density 
is independent to galaxy type and environment. We have found that core- like density profiles fit better the rotation 
curves than cuspy-like ones. No major differences have been found between field, cluster and compact group galaxies 
with respect to their dark halo density profiles. 

Based on a sample of « 100 velocity fields of spiral galaxies belonging to 25 HCGs, our previous study has demon- 
strated that: 40% (38 galaxies) of the HCG galaxies analysed in our sample do not allow the computation of the 
RC. This is because these galaxies are strongly perturbed by interaction and mergers. 33% (31 galaxies) of HCG 
galaxies are regular enough to calculate a RC but the RC is not symmetric to allow the derivation of a mass model. 
These perturbed galaxies are mildly interacting. The remaining 27% (25 galaxies) of them are suitable to derive mass 
models. For a variety reasons only 18 galaxies have been studied in the present work. The RCs have been combined 
with 2MASS J-band surface brightness profiles to derive their mass models. Two shapes of dark halos have been 
considered: a core- like density profile (isothermal sphere or ISO) and a cuspy-like one (NFW [4]). Best fit models and 
maximum disk models have been computed for the two dark halo profil! es. We have compared the HCG galaxies 
with two samples of galaxies in different environments: field and clusters galaxies. The three samples (HCG, isolated 
and cluster) have been analysed using the same tools. 

The results are summarized below: 

• No obvious differences can be found between the halo parameters for the HCG galaxies and galaxies in other 
environments. The strong correlation between po and ro is present for the three samples and the slopes of the 
linear regressions between them are very similar (between -1.62 and -1.00), well inside the uncertainties. In a 
study based on late-type and dwarf spheroidal galaxies, [5] confirms this correlation giving a slope of -1.038. On 
the other hand, the linear regression constant found by [5] is different in comparison with those found with the 
three samples. This scale difference is may be due to the fact that [5] compiled his sample using several sources 
using different methods to obtain halo parameters. 

• The use of the NFW model gives less satisfactory results than the ISO model. As mentioned by S08 [3], the 
coeflacient is usually larger when using the NFW. NFW model gives a worse correlation between the halo central 
density and the core radius than ISO. This is also true if we use the field galaxies sample. In contrast, the NFW 
model gives a consistent result for the cluster galaxies sample. Halo profiles are closer to isothermal spheres 
than NFW profile as already found by other authors [3], [6], [7]. The slope of the linear regression between pQ 
and ro is higher for the cluster galaxies sample than for the two other samples when the NFW model is used. 
The disagreement between ISO and NFW is smaller with cluster galaxies sample. 61% of HCG galaxies shows 
higher disk M/L using the ISO model compared to NFW. No relation between the disk scale length and the 
halo central density is seen using either ISO or NFW models. 

• We explored the possible connection between the halo parameters, the halo mass fraction and the disk M/L. 
The halo mass is high for both field galaxies and compact groups galaxies (75 to 95% of the total mass), leaving 
only modest room for the disk mass. The halo surface density is independent of the absolute B magnitude and 
no clear relation is seen between the disk M/L and Mb, meaning that the halo is independent of the galaxy 
luminosity. 

Recently, [8] suggest that core halos observed using high-resolution velocity fields in real dark matter-dominated 
galaxies are genuine and cannot be ascribed to systematic errors, halo triaxiality, or non-circular motions. 
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FIG. 17: Example of rotation curve and the fit of the different components: bulge, dislt and isothermal halo. 
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FIG. 18: Different correlation found with halos parameters and magnitude and disk scale length. 
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Warm dark matter: the linear Boltzmann-Vlasov equation vs. observations 

The distribution function of the decoupled dark matter (DM) particles can be written as 

where /(f (p) is the zeroth order freezed-out DM distribution function in or out of thermal equilibrium and Fi {x, p; t) 
the primordial fluctuations around it. The evolution of Fi is governed by the Boltzmann-Vlasov equation linearized 
around fo"^\p)- We evolve Fi{x,p; t) according to the linearized Boltzmann-Vlasov equation since the end of inflation 
till today. The DM density fluctuations are given by 

Ait, k)^m J ^ J d'x e-"-^ F,{x,p; t) . 

The initial data for the evolution are the primordial inflationary fluctuations 

Tig -1 

|0fe| = V2 7r^ ' , |Ao| ~ 4.94 10-^ ~ 0.964, fco = 2 Gpc-S 

A: 2 \Ko/ 

where g{k) is a random gaussian held. It is convenient to factorize the initial data as A(z, k) ~ pdm A(z, k) g(k) . 
We show that the linearized Boltzmann-Vlasov equation can be recasted as a pair of Volterra integral equations for 
Warm Dark Matter (WDM) plus the neutrino density fluctuations A{z,k) [1]. When the WDM becomes non- 
relativistic, the pair of integral equations reduces to the single equation: 

A{z,k) = h{z,k)+ ^ / ds'U{kru,M^)-s']}A{z{s'),k), (3.5) 

(Z + I) k run J so 

where z(s)+l = (zg^+l) sinh^ s , Zeq+l ~ 3200 , A(initial, k) = 1. Here h{z, k) is a known function: it contains 
the memory from previous ultra-relativistic (UR) evolution of the DM and n{x) = Q dQ f^^'{Q) sin(Q x) . 
The integral equation p.Sp is valid both in the radiation RD and matter MD dominated eras. Eq. (|3.5p becomes the 
Gilbert equation in the MD era (plus memory terms). 

The physical characteristic length scale in this linear evolution is the free streaming scale (or Jeans' scale) 

2_ 

( 3 J\/2 \ _ 1 X 

rur, = 2 y/l + JTq ^ fT^'^ = 21.1 qp ^ kpc , ~ 0.4 m/keV , 

\-tlQ ^/iljjM Qprirn J 

where qp = Qprim/O^^y)'^ is the primordial phase space density. DM particles can freely propagate over distances 
less or equal the free streaming scale. 

We plot in fig. [19] the transfer function A(z = 0, fc)/A(initial, k) vs. k run computed from ea. p.Sp for particles of 
mass TO = 1 keV. The curve in red corresponds to fermions decoupling at thermal equilibrium and the curve in blue 
to sterile neutrinos in the x-model decoupling out of thermal equilibrium. 

The linear density profile pun{r, z) at redshift z is given by 

PUn{r,z) ^ ——— I k dk sm{kr) A{k,z) for g(fc) = 1 
2 7r^ r Jq 

Notice that A{z,j) ~ A(0,7). Therefore, the profile shape turns to be redshift independent in the MD/A era. 
We plot in fig. [501 the linear density profiles piin{f) / Pun{^) vs. x = r/run- WDM linear density profiles turn to be 
cored for r — >■ 0. We find in the intermediate regime r > run'. 

Pi,„(r) ~= Co (^j PH„(0) , H-n,/2 = 1.482 

Namely, pun (r) scales with the primordial spectral index Us . These theoretical linear results agree with the universal 
empirical behaviour 7'-i-6±o.4. q Walker et al. (2009) (observations) and with I. M. Vass et al. (2009) (sim- 
ulations). In the asymptotic regime r 3> Hin, the small k behaviour of A(fc,itoday) ci (fc runY with s ~ 0.5 
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FIG. 19: The WDM transfer function A(z = 0, A:)/A(initial, A:) vs. k run computed from ea. (|3.5|l for 1 keV particles. Red: 
tliermal equilibrium decoupling. Blue: out of equilibrium decoupling. 




FIG. 20: Linear density profiles punij') / piin{Q) vs. x = r/run. Same colur code as in fig. [191 

implies the presence of a tail: piinij) — c (run/r) . For WDM, the agreement between the linear theory and 
the observations is remarkable. Linear CDM profiles turn to be cusped: PHn(j')cDM vs. r exhibits a cusp behaviour 
for r > 0.03 pc, (in agreement with CDM simulations, which show the predictivity of our results). Notice that 
observations in DM dominated galaxies always find cores as predicted by WDM. 

The expected overdensity within a radius R at redshift z in the linear regime is given by 

a^{R, 2) = / — A2(z, k) W^{kR) , W{kR) = window function . 
Jo k 

(J^{R, z) ~ 1 is the borderline between linear and non-linear regimes (sec fig. [2T|) . Objects (galaxies) of scale R and 
mass ^ R^ start to form when this scale becomes non-linear. Our linear claculation shows that smaller objects form 
earlier in agreement with astronomical observations [2] . 

In short, WDM is characterized by : (i) its initial power spectrum cutted off for scales below run ^ 50 kpc as shown 
by Fig. [m Thus, structures are not formed in WDM for scales below run ^ 50 kpc. (ii) its initial velocity dispersion. 
However, this is negligible for z < 20 where the non-linear regime starts. 
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FIG. 21: log M/Mq vs. \og{z + 1) for CDM (violet) and WDM: 1 keV (red), 2 keV (green), 4 keV (blue) DM particles 
decoupling in equilibrium and 1 keV (light-blue) sterile neutrinos 
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Perturbation theories for large-scale structures 

In the standard cosmological scenario, the large-scale structures we observe in the recent Universe (galaxies, clusters, 
filaments, voids, ..) have formed through the amplification by gravitational instability of small primordial perturba- 
tions generated in the primordial Universe by quantum fluctuations. Since the initial power increases on small scales 
(for popular models such as CDM), on large scales or at early times it is possible to use linear theory or, more generally, 
perturbation theory, while on small scales or at late times one must use numerical simulations or phcnomenological 
models such as the halo model. 

In this talk, we have described some properties of perturbative approaches, in particular the different behaviors found 
in Eulerian and Lagrangian frameworks, and we have shown how to combine these schemes with phcnomenological 
halo models, to build unified models that can describe all scales. 

Focusing on the power spectrum P{k), we have first recalled that standard perturbation theory is not well-behaved, 
because higher-order terms grow increasingly fast at high k with strong cancellations between various orders. On the 
other hand, it can be checked that many orders of perturbation theory (typically up to order 30 at 2: = 2) are relevant 
before non-pcrturbativc corrections associated with shell crossing dominate [1]. Therefore, it is useful to build other 
perturbative schemes (i.e. partial rcsummations) that are better behaved and more efficient. 




k [h Mpc-iJ k [hMpc-iJ 



FIG. 22: Left panel: ratio of the nonlinear power spectrum P{k) to a smooth linear power spectrum Pls without acoustic 
baryonic oscillation, at z = 1. The points with error bars are the results from N-body simulations, the black solid line is 
the linear power spectrum, the upper blue dashed line is the prediction of standard perturbation theory up to 1-loop order, 
the green solid line is the theoretical prediction (combining the "large-N" perturbative resummation and the halo model), the 
magenta dot-dashed line is a popular fit to simulations [9]. Right panel: bispectrum for equilateral configurations, at z — 1. 
The points are the results from numerical simulations, the black solid line is the tree-order result, the blue dot-dashed line is 
the standard 1-loop order prediction, the red solid line is the full theoretical prediction that also includes the non-perturbative 
"1-halo" contribution. The magenta dot-dashed line is a simple phenomenological model [10]. 

As a first example, wc have described a "large-N" perturbative approach, based on a path-integral formalism. At 
one-loop order this provides a significantly higher accuracy on large scales than the standard perturbation theory, 
while being well-behaved at high k [2]. Next, we have explained how an alternative approach [3], based on a high-fc 
resummation and the assumption of a wide separation of scales, corresponds to replacing the equations of motion by 
linear equations with random coefficients. This yields Eulerian-space propagators that show a Gaussian decay at high 
k, in good agreement with simulations. However, this behavior is not a true loss of memory of the density field, but 
a "sweeping effect" due to long wavelength modes of the velocity field that coherently move density structures [4] . 
In particular, applying the same approach in Lagrangian-space, which is not sensitive to this "sweeping effect" , one 
obtains a propagator that no longer decays [5] . 

Then, we have shown how these two distinct behaviors can be understood from the study of the "adhesion model" , 
where exact results can be derived. This model, introduced by [6], extends the Zcldovich dynamics by making particles 
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glue together at collisions. This dynamics builds a cosmic web and density fields that are very similar to those found 
in cosmological simulations, except that halos are pointlike masses. Then, one can show that the Eulerian-space 
propagator is actually equal to the velocity probability distribution, with an exponential-like high-fc decay, whereas 
the Lagrangian propagator can be written in terms of the halo mass function, with a power- law high-fc decay [7]. This 
confirms the explanation of these two different behaviors and explicitly shows that Lagrangian propagators should be 
more sensitive probes of the density field than their Eulerian counterparts. 

Next, we have described how to build unified models than combine perturbation theories with halo models. More 
precisely, focusing on the power spectrum P{k), the "2-halo" contribution can be written in terms of the perturbative 
prediction, whereas the non-perturbative "1-halo" term involves the halo mass function and density profile. Moreover, 
we have pointed out that the latter contains a specific countcrtcrm that ensures a physically meaningful behavior on 
large scales and is needed to reach a high accuracy [2] . 

Then, we have described how the same approach allows us to obtain the density bispectrum, from large to small 
scales [8]. 



Finally, we have briefly reviewed a few other perturbative approaches and discussed the Vlasov-Poisson system, 
which may allow going beyond the fluid approximation and taking into account some shell-crossing effects. 
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Supernova bound on keV-mass sterile neutrinos 

Sterile neutrinos with masses in the keV range, which are a promising candidate for warm dark matter, can be 
copiously produced in the supernova (SN) core. For m,, > 100 keV, the vacuum mixing angle of sterile neutrinos is 
stringently constrained sin^ 29 < 10"^ in order to avoid excessive energy loss [1-5]. For smaller masses, however, the 
Mikheyev-Smirnov-Wolfenstein (MSW) matter effect on active-sterile neutrino mixing becomes very important and 
the SN bound on vacuum mixing angle is not that obvious. Note that the bounds on mixing angles depend on which 
flavor the sterile neutrino mixes with. We concentrate on the SN bound in i/T--j/s-mixing case for simplicity, because 
Ut and Ft- only have neutral-current interactions and essentially stay in thermal equilibrium with ambient matter. 

The matter density in the SN core is so high that the incoherent scattering of active neutrinos on matter particles 
may even dominate over flavor oscillations as the production mechanism for keV-mass sterile neutrinos. An elegant 
formalism to deal with both incoherent scattering and flavor oscillations have been developed in Refs. [6-8], where the 
evolution equations for the occupation numbers of different neutrino species have been derived. In the weak-damping 
limit, which is always valid for supernova neutrinos mixing with keV-mass sterile neutrinos, the evolution of i^r number 
density is determined by [9] 

where S2e„ = sin20y with being the neutrino mixing angle in matter, the occupation number of and W^,^ 
the transition probability for i'{E') -t- a — >■ ^{E) + a with a being background particles in the SN core. In a similar 
way, we can derive the evolution equation of the P,- number density, involving the mixing angle 9p, the occupation 
number and the transition probability W^,^. Due to the MSW effect, the mixing angle of neutrinos in matter is 
different from that of antineutrinos 

sin2 20, p = ""-^^ , (3.7) 

' &in^2e + {co&2e±E/E,Y ' ^ ' 

where 6 denotes the vacuum mixing angle, and the upper sign refers to v and the lower to v. The resonant energy 
E, = AmV2|Kj is defined as 



E, = 3.25 MeV [-^^) Pil \Yo - Y^X^ , (3.8) 



where pi4 is the matter density p in units of 10^^ g cm^'^ and Yq = (1 — Yg — '^Y,^)/4:. Note that Yx = {Nx — Nx)/Nb 
with TVb being the baryon number density, Nx and Nx being the number densities of particle x and its antiparticle x. 
For tan neutrinos, the matter potential V,^ — ~{Gf/^/2)Nb (1 — Y^ — 2Y,^ — 4Y,^) is negative if the typical values 
of Ye = 0.3, Yij^ ~ 0.07 and Y,^ = for a SN core are taken. Therefore, the mixing angle for i>r is enhanced by 
matter effects, and the emission rate for D^. exceeds that for 1/^, indicating that a Vt-Vt asymmetry (i.e., Y,^^ ^ 0) 
will be established. An interesting feedback effect emerges: (i) The chemical potential for tau neutrinos develops 
and thus changes the occupation numbers of and Vt] (ii) The Vt-v-t asymmetry shifts the resonant energy iS^, 
and thus modifies the mixing angles dy and (iii) Both effects in (i) and (ii) will feed back on the emission rates. 
Hence a stationary state of this active-sterile neutrino system will be achieved if the emission rates for neutrinos and 
antineutrinos become equal to each other [9]. 

Given the sterile neutrino mass rris and vacuum mixing angle 0, the energy loss rate E(t) due to sterile neutrino 
emission can be calculated by following the evolution of v^-v^ asymmetry Y,^ (i). It has been found that the stationary 
state can be reached within one second and the feedback effect is very important for 20 keV < < 80 keV and 
10~^ < sin^ 2d < 10~^. To avoid excessive energy losses, we require that the average energy-loss rate {£) = J^"^ £{t) di 
with Td = 1 s should be {£) < 3.0 x 10'^'' erg cm~^ s~^. Otherwise, the duration of neutrino burst from Supernova 
1987A would have been significantly reduced. In Fig. [23l we show the contours of energy- loss rates in the (sin^ 29, ms)- 
plane, where we have assumed a homogeneous and isotropic core with matter density p = 3.0 x 10~^^ g cm~'^ and 
temperature T = 30 MeV. Based on the energy-loss argument, the purple region has been excluded. The most 
stringent bound sin^ 29 < 10~^ arises for rus = 50 keV. For the large-mixing angle region, the energy-loss rate is 
actually small, because sterile neutrinos have been trapped in the core and cannot carry away energies. However, 
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FIG. 23: Supernova bound on sterile neutrino masses ttis and mixing angles 6, where the purple region is excluded by the 
energy-loss argument while the green one by the energy-transfer argument [9]. The excluded region will be extended to the 
dashed (red) line if the build-up of degeneracy parameter is ignored, i.e., ri(t) = 0. The dot-dashed (green) line represents the 
sterile neutrinos as dark matter with the correct relic abundance Q,sh? — 0.1. 



the mean free path of sterile neutrinos is comparable to or even longer than that of ordinary neutrinos, indicating 
that they may transfer energies in a more efficient way. As a consequence, the duration of neutrino burst will be 
shortened by emitting neutrinos more rapidly. In this sense, the excessive energy transfer should be as dangerous as 
the excessive energy loss. Hence the large-mixing angle region is indeed excluded if the energy-transfer argument is 
applied. The green line in Fig. indicates the relic abundance of dark matter Q.sh'^ = 0.1, where keV-mass sterile 
neutrinos are warm dark matter and the non-resonant production mechanism is assumed. If we ignore the feedback 
effect (i.e., a vanishing chemical potential for tau neutrinos rj — /i^^/T = 0), the excluded region will extend to the 
red line, which overlaps the relic-abundance line. As shown in Fig. I23[ however, the mixing angles are essentially 
unconstrained in the favored warm-dark-matter mass range 1 keV < < 10 keV [10,11]. 

As for the i/^-i/j- mixing case, our discussions about the feedback effects are essentially applicable. However, the 
charged-current interactions of and should be taken into account, and the change of I'^-i^^ asymmetry will be 
redistributed between muon neutrinos and charged muons. The Ve-i^s mixing in SN cores is more involved because of 
the large trapped electron number and high j/g degeneracy. Besides energy loss, deleptonization by sterile neutrino 
emission is an effect to be taken into account. It has been suggested that keV-mass sterile neutrinos mixing with 
electron neutrinos may help supernova explosions [12,13]. However, this case requires a dedicated investigation. 
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IV. SUMMARY AND CONCLUSIONS OF THE WORKSHOP BY H. J. DE VEGA AND 

N. G. SANCHEZ 

A. General view and clarifying remarks 

Participants came from Em'ope, North and South America, Russia, Japan, Korea. Discussions and lectures were 
outstanding. WDM research evolves fastly in both astronomical, numerical, theoretical, particle and experimental 
research. The Workshop allowed to make visible the work on WDM made by different groups over the world and 
cristalize WDM as the viable component of the standard cosmological model in agreement with CMB + Large Scale 
Structure (LSS) + Small Scale Structure (SSS) observations, AWDM, in contrast to ACDM which only agree with 
CMB+LSS observations and is plagued with SSS problems. 

The participants and the programme represented the different communities doing research on dark matter: 

• Observational astronomers 

• Computer simulators 

• Theoretical astrophysicists not doing simulations 

• Physical theorists 

• Particle experimentalists 

WDM refers to kcV scale DM particles. This is not Hot DM (HDM). (HDM refers to cV scale DM particles, which 
are already ruled out). CDM refers to heavy DM particles (so called wimps of GcV scale or any scale larger than 
keV). 

It should be recalled that the connection between small scale structure features and the mass of the DM particle 
follows mainly from the value of the free-streaming length / fg. Structures smaller than / fg are erased by free-streaming. 
WDM particles with mass in the keV scale produce Ifs ^ 100 kpc while 100 GeV CDM particles produce an extremely 
small Ifs ^ 0.1 pc. While the keV WDM Ifs ^ 100 kpc is in nice agreement with the astronomical observations, the 
GeV CDM Ifs is a million times smaller and produces the existence of too many small scale structures till distances 
of the size of the Oort's cloud in the solar system. No structures of such type have ever been observed. 

Also, the name CDM precisely refers to simulations with heavy DM particles in the GcV scale. Most of the literature 
on CDM simulations do not make explicit the relevant ingredient which is the mass of the DM particle (GeV scale 
wimps in the CDM case). 

The mass of the DM particle with the free-streaming length naturally enters in the initial power spectrum used in 
the N-body simulations and in the initial velocity. The power spectrum for large scales beyond 100 kpc is identical 
for WDM and CDM particles, while the WDM spectrum is naturally cut off at scales below 100 kpc, corresponding 
to the keV particle mass free-streaming length. In contrast, the CDM spectrum smoothly continues for smaller and 
smaller scales till ~ 0.1 pc, which gives rise to the overabundance of predicted CDM structures at such scales. 

CDM particles are always non-relativistic, the initial velocities are taken zero in CDM simulations, (and phase 
space density is unrealistically infinity in CDM simulations), while all this is not so for WDM. 

Since keV scale DM particles are non rclativistic for z < 10^ they could also deserve the name of cold dark matter, 
although for historical reasons the name WDM is used. Overall, seen in perspective today, the reasons why CDM 
does not work are simple: the heavy wimps are excessively non-relativistic (too heavy, too cold, too slow), and thus 
frozen, which preclude them to erase the structures below the kpc scale, while the eV particles (HDM) are excessively 
rclativistic, too light and fast, (its free streaming length is too large), which erase all structures below the Mpc scale; 
in between, WDM keV particles produce the right answer. 



B. Conclusions 



Some conclusions are: 
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Sterile neutrinos with mass in the keV scale (1 to 10 keV) emerge as leading candidates for the dark matter 
(DM) particle from theory combined with astronomical observations. 

DM particles in the keV scale (warm dark matter, WDM) naturally reproduce (i) the observed galaxy structures 
at small scales (less than 50 kpc), (ii) the observed value of the galaxy surface density and phase space density 
(iii) the cored profiles of galaxy density profiles seen in astronomical observations. 

Heavier DM particles (as wimps in the GeV mass scale) do not reproduce the above important galaxy obser- 
vations and run into growing and growing serious problems (they produce satellites problem, voids problem, 
galaxy size problem, unobserved density cusps and other problems). 

Sterile neutrinos are serious WDM candidates: Minimal extensions of the Standard Model of particle physics 
include keV sterile neutrinos which are very weakly coupled to the standard model particles and are produced 
via the oscillation of the light (eV) active neutrinos, with their mixing angle governing the amount of generated 
WDM. The mixing angle theta between active and sterile neutrinos should be in the lO""* scale to reproduce 
the average DM density in the Universe. 

Sterile neutrinos are usually produced out of thermal equilibrium. The production can be non-resonant (in the 
absence of lepton asymmetries) or resonantly ennhanced (if lepton asymmetries are present). The usual X ray 
bound together with the Lyman alpha bound forbids the non-resonant mechanism in the i^MSM model. 

Sterile neutrinos can decay into an active-like neutrino and an X-ray photon. Abundance and phase space 
density of dwarf spheroidal galaxies constrain the mass to be in the ~ keV range. Small scale aspects of sterile 
neutrinos and different mechanisms of their production were presented: The transfer function and power spectra 
obtained by solving the coUisionlcss Boltzmann equation during the radiation and matter dominated eras feature 
new WDM acoustic oscillations on mass scales ^ 10® — 10^ Mq. 

Lyman-alpha constraints have been often misinterpreted or superficially invoked in the past to wrongly suggest 
a tension with WDM, but those constraints have been by now clarified and relaxed, and such a tension does not 
exist: keV sterile neutrino dark matter (WDM) is consistent with Lyman-alpha constraints within a wide range 
of the sterile neutrino model parameters. Only for sterile neutrinos assuming a non-resonant (Dodclson- 
Widrow model) production mechanism, Lyman-alpha constraints provide a lower bound for the mass of about 
4 keV . For thermal WDM relics (WDM particles decoupling at thermal equilibrium) the Lyman-alpha lower 
particle mass bounds are smaller than for non-thermal WDM relics (WDM particles decoupling out of thermal 
equilibrium). The number of Milky- Way satellites indicates lower bounds between 1 and 13 keV for different 
models of sterile neutrinos. 

WDM kcV sterile neutrinos can be copiously produced in the supernovae cores. Supcrnovae stringently constraint 
the neutrino mixing angle squared to be < 10~^ for sterile neutrino masses m > 100 keV (in order to avoid 
excessive energy lost) but for smaller sterile neutrino masses the SN bound is not so direct. Within the models 
worked out till now, mixing angles are essentially unconstrained by SN in the favoured WDM mass range, 
namely 1 < m < 10 keV. Mixing between electron and keV sterile neutrinos could help SN explosions, case 
which deserve investigation 

Signatures for a right-handed few keV sterile neutrino should be Lyman alpha emission and absorption at around 
a few microns; corresponding emission and absorption lines might be visible from molecular Hydrogen H2 and H3 
and their ions in the far infrared and sub-mm wavelength range. The detection at very high redshift of massive 
star formation, stellar evolution and the formation of the first super-massive black holes would constitute the 
most striking and testable prediction of WDM sterile neutrinos. 

The effect of kcV WDM can be also observable in the statistical properties of cosmological Large Scale Structure. 
Cosmic shear (weak gravitational lensing) does not strongly depend on baryonic physics and is a promisingig 
probe. First results in a simple thermal relic scenario indicate that future weak lensing surveys could see a 
WDM signal for itlwdm ^ 2 kcV or smaller. The predicted limit beyond which these surveys will not see a 
WDM signal is mwDM ~ 2.5 kcV (thermal relic) for combined Euclid -I- Planck. More realistic models deserve 
investigation and are expected to relaxe such minimal bound. With the real data, the non-linear WDM model 
should be taken into account 

The possibility of laboratory detection of warm dark matter is extremely interesting. Only a direct detection 
of the DM particle can give a clear-cut answer to the nature of DM. At present, only the Katrin and Mare 
experiments have the possibility to do that for sterile neutrinos. Mare bounds on sterile neutrinos have been 
reported in this Workshop. Namely, bounds from the beta decay of Rel87 and EC decay of Hol63. Mare keeps 
collecting data in both. 
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The possibility that Katrin experiment can look to sterile neutrinos in the tritium decay did appeared in the 
discussions. Katrin experiment have the potentiality to detect warm dark matter if its set-up would be adapted 
to look to keV scale sterile neutrinos. Katrin experiment concentrates its attention right now on the electron 
spectrum near its end-point since its goal is to measure the active neutrino mass. Sterile neutrinos in the tritium 
decay will affect the electron kinematics at an energy about m below the end-point of the spectrum (m = sterile 
neutrinos mass). Katrin in the future could perhaps adapt its set-up to look to keV scale sterile neutrinos. It 
will be a a fantastic discovery to detect dark matter in a beta decay 

Astronomical observations strongly indicate that dark matter halos are cored till scales below 1 kpc. 

More precisely, the measured cores are not hidden cusps. CDM Numerical simulations -with wimps (particles 
heavier than 1 GeV)- without and with baryons yield cusped dark matter halos. Adding baryons do not alleviate 
the problems of wimps (CDM) simulations, on the contrary adiabatic contraction increases the central density 
of cups worsening the discrepancies with astronomical observations. In order to transform the CDM cusps into 
cores, the baryon-l-CDM simulations need to introduce strong baryon and supernovae feedback which produces 
a large star formation rate contradicting the observations. None of the predictions of CDM simulations at 
small scales (cusps, substructures, ...) have been observed. The discrepancies of CDM simulations with the 
astronomical observations at small scales < 100 kpc is staggering: satellite problem (for example, only 1/3 
of satellites predicted by CDM simulations around our galaxy are observed), the surface density problem (the 
value obtained in CDM simulations is 1000 times larger than the observed galaxy surface density value), the 
voids problem, size problem (CDM simulations produce too small galaxies). 

The use of keV scale WDM particles in the simulations instead of the GeV CDM wimps, alleviate all the above 
problems. For the core-cusp problem, setting the velocity dispersion of keV scale DM particles seems beyond the 
present resolution of computer simulations. However, the velocity dispersion of WDM particles is negligible for 
2; < 20 where the non-linear regime and structure formation starts. Analytic work in the linear approximation 
produces cored profiles for keV scale DM particles and cusped profiles for CDM. Model-independent analysis of 
DM from phase-space density and surface density observational data plus theoretical analysis points to a DM 
particle mass in the keV scale. The dark matter particle candidates with high mass (100 GeV, 'wimps') are 
strongly disfavored, while cored (non cusped) dark matter halos and warm (keV scale mass) dark matter are 
strongly favoured from theory and astrophysical observations. 

An 'Universal Rotation Curve' (URC) of spiral galaxies emerged from 3200 individual observed Rotation Curves 
(RCs) and reproduces remarkably well out to the virial radius the Rotation Curve of any spiral galaxy. The 
URC is the observational counterpart of the circular velocity profile from cosmological simulations. ACDM 
simulations give the well known NFW cuspy halo profile. A careful analysis from about 100 observed high 
quality rotation curves has now ruled out the disk -I- NFW halo mass model, in favor of cored profiles. 
The observed galaxy surface density (surface gravity acceleration) appears to be universal within ^ 10% with 
values around 120 Mq/pc^ , irrespective of galaxy morphology, luminosity and Hubble types, spanning over 14 
magnitudes in luminosity and mass profiles determined by several independent methods. 

Interestingly enough, a constant surface density (in this case column density) with value around 120 Mq/pc^ 
similar to that found for galaxy systems is found too for the interstellar molecular clouds, irrespective of size and 
compositions over six order of magnitude; this universal surface density in molecular clouds is a consequence of 
the Larson scaling laws. This suggests the role of gravity on matter (whatever DM or baryonic) as a dominant 
underlying mechanism to produce such universal surface density in galaxies and molecular clouds. Recent re- 
examination of different and independent (mostly millimeter) molecular cloud data sets show that intcrcstcUar 
clouds do follow Larson law Mass ^ (Size)'^ exquisitely well, and therefore very similar projected mass densities 
at each extinction threshold. Such scaling and universality should play a key role in cloud structure formation. 

Visible and dark matter observed distributions (from rotation curves of high resolution 2D velocity fields) of 
galaxies in compact groups, and the comparison with those of galaxies in clusters and field galaxies show that: (i) 
The central halo surface density is constant with respect to the total absolute magnitude similar to what is found 
for the isolated galaxies, suggesting that the halo density is independent of the galaxy type and environment, 
(ii) Core density profiles fit better the rotation curves than cuspy profiles. Dark halo density profiles are found 
almost the same in field galaxies, cluster and compact group galaxies. Core halos observed using high-resolution 
velocity fields in dark matter galaxies are genuine and cannot be ascribed to systematic errors, halo triaxiality, 
or non-circular motions, (iii) The halo mass is high (75 to 95 % of the total mass) for both field galaxies and 
compact groups galaxies, living modest room for a dark mass disk. No relation between the disk scale length 
and the halo central density is seen, the halo being independent of galaxy luminosity. 
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Besides the DM haloes, ACDM models ubiquituosly predict eold dark matter disks, [also believed as the conse- 
quence of merging events]: CDM simulations predict that a galaxy such as the Milky Way should host a CDM 
disk (with a scale height of 2.1 — 2.4 kpc and a local density at the solar position 0.25 — 1.0 times that of the DM 
halo; some models reach a height of 4 kpc from the galactic plane), which is thicker than any visible disk (the 
scale height of the galactic old thick stellar disk is ^ 0.9 kpc; young stars and ISM form even thinner structures 
of height 0.3 and 0.1 kpc respectively). However, careful astronomical observations performed to see such disk 
have found no evidence of such CDM disk in the Milky Way. 

A key part of any galaxy formation process and evolution involves dark matter. Cold gas accretion and mergers 
became important ingredients of the CDM models but they have little observational evidence. DM properties 
and its correlation with stellar masses are measured today up to z = 2; at z > 2 observations are much less 
certain. Using kinematics and star formation rates, all types of masses -gaseous, stellar and dark- are measured 
now up to z = 1.4. The DM density within galaxies declines at higher redshifts. Star formation is observed to 
be more common in the past than today. More passive galaxies are in more massive DM halos, namely most 
massive DM halos have lowest fraction of stellar mass. CDM predicts high overabundance of structure today 
and under-abundance of structure in the past with respect to observations. The size-luminosity scaling relation 
is the tightest of all purely photometric correlations used to characterize galaxies; its environmental dependence 
have been highly debated but recent findings show that the size-luminosity relation of nearby elliptical galaxies is 
well defined by a fundamental line and is environmental independent. Observed structural properties of elliptical 
galaxies appear simple and with no environmental dependence, showing that their growth via important mergers 
-as required by CDM galaxy formation- is not plausible. Moreover, observations in brightest cluster galaxies 
(BCGs) show little changes in the sizes of most massive galaxies since z = 1 and this scale-size evolution appears 
closer to that of radio galaxies over a similar epoch. This lack of size growth evolution, a lack of BCG stellar 
mass evolution is observed too, demonstrates that major merging is not an important process. Again, 
these observations put in serious trouble CDM 'semianalytical models' of BCG evolution which require about 
70% of the final BCG stellar mass to be accreted in the evolution and important growth factors in size of massive 
elliptical massive galaxies. 

Recent AWDM N-body simulations have been performed by different groups. High resolution simulations for 
different types of DM (HDM, WDM or CDM), allow to visualize the effects of the mass of the corresponding DM 
particles: free-streaming length scale, initial velocities and associated phase space density properties: for masses 
in the eV scale (HDM), halo formation occurs top down on all scales with the most massive haloes collapsing 
first; if primordial velocities are large enough, free streaming erases all perturbations and haloes cannot form 
(HDM) . The concentration- mass halo relation for mass of hundreds eV is reversed with respect to that found for 
CDM wimps of GeV mass. For realistic keV WDM these simulations deserve investigation: it could be expected 
from these HDM and CDM effects that combined free-streaming and velocity effects in keV WDM simulations 
could produce a bottom-up hierarchical scenario with the right amount of sub-structures (and some scale at 
which transition from top-down to bottom up regime is visualized). 

Moreover, interestingly enough, recent large high resolution AWDM N-body simulations allow to discriminate 
among thermal and non-thermal WDM (sterile neutrinos): Unlike conventional thermal relics, non-thermal 
WDM has a peculiar velocity distribution (a little skewed to low velocities) which translates into a character- 
istic linear matter power spectrum decreasing slowler across the cut-off free-streaming scale than the thermal 
WDM spectrum. As a consequence, the radial distribution of the subhalos predicted by WDM sterile neutrinos 
remarkably reproduces the observed distribution of Milky Way satellites in the range above ~ 40 kpc, while the 
thermal WDM supresses subgalactic structures perhaps too much, by a factor 2 — 4 than the observation. Both 
simulations were performed for a mass equal to 1 keV. Simulations for a mass larger than 1 keV (in the range 
between 2 and 10 keV, say) should still improve these results. 

The predicted AWDM galaxy distribution in the local universe (as performed by CLUES simulations with a 
a mass of mwDM = 1 keV) agrees well with the observed one in the ALFALFA survey. On the contrary, 
ACDM predicts a steep rise in the velocity function towards low velocities and thus forecasts much more sources 
than the ones observed by the ALFALFA survey (both in Virgo-direction as well as in anti- Virgo-direction) . 
These results show again the ACDM problems, also shown in the spectrum of mini-voids. AWDM provides a 
natural solution to these problems. WDM physics effectively acts as a truncation of the ACDM power spectrum. 
AWDM CLUES simulations with 1 kcV particles gives much better answer than ACDM when reproducing sizes 
of local minivoids. The velocity function of 1 keV WDM Local Volume-counterpart reproduces the observational 
velocity function remarkably well. Overall, keV WDM particles deserve dedicated experimental detection efforts 
and simulations. 
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The features observed in the cosmic-ray spectrum by Auger, Pamela, Fermi, HESS, CREAM and others can 
be all quantitatively well explained with the action of cosmic rays accelerated in the magnetic winds of very 
massive star explosions such as Wolf-Rayet stars, without any significant free parameter. All these observations 
of cosmic ray positrons and electrons and the like are due to normal astrophysical sources and processes, and do 
not require an hypothetical decay or annihilation of heavy CDM particles (wimps). The models of annihilation 
or decay of heavy CDM wimps are highly tailored to explain these normal astrophysical processes and their 
ability to survive observations is more than reduced. 

Theoretical analytic perturbative approachs for large scales in Eulerian and Lagrangian frameworks can be com- 
bined with phenomenological halo models to build unified schemes for describe all scales. The large k Gaussian 
decay of Eulerian-space propagators (in agreement with simulations) is not a true loss of memory of the density 
field but a "sweepping" effect of the velocity field modes which coherently move density structures. Lagrangian- 
space propagators are not sensitive to such effect. Extending the Zeldovich approximation ("adhesion model") 
show that Eulerian propagators are in fact the velocity probability distribution, whereas Lagrangian propagators 
explicitely relate to the halo mass function and are most sensitive probes of the density field than their Eulerian 
conterparts. Extensions to the Vlasov-Poisson system deserve investigation. 

Theoretical analytic modelling of the halo mass function on small scales gains insight with the use of the 
mathematical excursion set theory in a path integral formulation, which allows direct comparison with numer- 
ical simulations and observational results. Stochastic modeling of the halo collapse conditions can be easily 
implemented in this formalism. 

As an overall conclusion, CDM represents the past and WDM represents the future in the DM research. CDM 
research is 20 years old. CDM simulations and their proposed baryonic solutions, and the CDM wimp candidates 
(~ 100 GeV) are strongly pointed out by the galaxy observations as the wrong solution to DM. Theoretically, 
and placed in perspective after more than 20 years, the reason why CDM does not work appears simple and clear 
to understand and directly linked to the excessively heavy and slow CDM wimp, which determines an excessively 
small (for astrophysical structures) free streaming length, and unrealistic overabundance of structures at these 
scales. On the contrary, new keV WDM research, keV WDM simulations, and keV scale mass WDM particles 
are strongly favoured by galaxy observations and theoretical analysis, they naturally work and agree with the 
astrophysical observations at all scales, (galactic as well as cosmological scales). Theoretically, the reason why 
WDM works so well is clear and simple, directly linked to the keV scale mass and velocities of the WDM 
particles, and free-streaming length. The experimental search for serious WDM particle candidates (sterile 
neutrinos) appears urgent and important: it will be a fantastic discovery to detect dark matter in a beta decay. 
There is a formidable WDM work to perform ahead of us, these highlights point some of the directions where 
it is worthwhile to put the effort. 



C. The present context and future in DM research. 

Facts and status of DM research: Astrophysical observations point to the existence of DM. Despite of that, 
proposals to replace DM by modifing the laws of physics did appeared, however notice that modifying gravity 
spoils the standard model of cosmology and particle physics not providing an alternative. After more than twenty 
active years the subject of DM is mature, (many people is involved in this problem, different groups perform 
N-body cosmological simulations and on the other hand direct experimental particle searches are performed 
by different groups, an important number of conferences on DM and related subjects is held regularly). DM 
research appears mainly in three sets: (a) Particle physics DM model building beyond the standard model of 
particle physics, dedicated laboratory experiments, annhilating DM. all concentrated on CDM and CDM wimps, 
(b) Astrophysical DM: astronomical observations, astrophysical models, (c) Numerical CDM simulations. : The 
results of (a) and (b) do not agree and (b) and (c) do not agree neither at small scales. None of the small scale 
predictions of CDM simulations have been observed: cusps and over abundance of substructures differ by a 
huge factor with respect to those observed. In addition, all direct dedicated searchs of CDM wimps from more 
than twenty years gave null results. Something is going wrong in the CDM research and the right answer is: the 
nature of DM is not cold (GeV scale) but warm (keV scale). 

Many researchers continue to work with heavy CDM candidates (mass > 1 GeV) despite the staggering 
evidence that these CDM particles do not reproduce the small scale astronomical observations (< 100 kpc). 
Why? [It is known now that the keV scale DM particles naturally produce the observed small scale structure] . 
Such strategic question is present in many discussions, everyday and off of the record (and on the record) talks in 
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the field. The answer deals in large part with the inertia (material, intellectual, social, other, ...) that structured 
research and big-sized science in general do have, which involve huge number of people, huge budgets, longtime 
planned experiments, and the "power" (and the conservation of power) such situation could allow to some of the 
research lines following the trend; as long as budgets will allow to run wimp experimental searches and CDM 
simulations such research lines could not deeply change, although they would progressively decline. 

• Notice that in most of the DM litterature or conferences, wimps are still "granted" as "the" DM particle, and 
CDM as "the" DM; is only recently that the differences and clarifications are being clearly recognized and 
acknowledged. While wimps were a testable hypothesis at the beginning of the CDM research, one could ask 
oneself why they continue to be worked out and "searched" experimentally in spite of the strong astronomical 
and astrophysical evidence against them. 

Similar situations (although not as extremal as the CDM situation) happened in other branches of physics and 
cosmology: Before the CMB anisotropy observations, the issue of structure formation was plugged with several 
alternative proposals which were afterwards ruled out. Also, string theory passed from being considered "the 
theory of everything" to "the theory of nothing" (as a physical theory), as no physical experimental evidence 
have been obtained and its cosmological implementation and predictions desagree with observations. (Despite all 
that, papers on such proposals continue -and probably will continue- to appear. But is clear that big dedicated 
experiments are not planned or built to test such papers). In science, what is today 'popular' can be discarded 
afterwards; what is today 'new' and minoritary can becomes 'standard' and majoritarily accepted if verified 
experimentally. 



'Examine the objects as they are and you will see their true nature; look at them from your own ego 
and you will see only your feelings; because nature is neutral, while your feelings are only prejudice 

and obscurity' 

[Gerry Gilmore quoting Shao Yong, 1011-1077 in the 14th Paris Cosmology Colloquium Chalonge 2010 



http: / / chalonge.obspm.fr /Programme JF'aris2010. html larXiv: 1 009 . 3494] . 



The Lectures of the Workshop can be found at: 



http:/ /www. chalonge. obspm.fr/Programme_CIAS2011. html 

The photos of the Workshop can be foimd at: 



http: / /www. chalonge. obspm.fr/albumCIAS2011/index.html 
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